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Abstract. We have compiled a sample of disk galaxies with available photometry in the B and K bands, velocity 
line- widths and HI integral fluxes. Several parameters which trace the luminous, baryonic and dark matter contents 
were inferred. We investigated how these parameters do vary with different galaxy properties, and confronted 
the results with predictions of galaxy evolutionary models in the context of the A Cold Dark Matter (ACDM) 
cosmogony. The ratio of disk-to-total maximum circular velocity, (Vd, m /Vt,m), depends mainly on the central 
disk surface density Ed,o (or surface brightness, SB), increasing roughly as Ed,o 0,15 - While a fraction of high SB 
galaxies have a (Vd, m /K,m) ratio corresponding to the maximum disk solution, the low SB are completely dark 
matter dominated. The trend is similar for the models, although they have slightly smaller (Vd, m /K,m) ratios 
than observations, in particular at the highest SBs and when small baryon fractions are used. The scatter in the 
(Vd,m/Vt,m) — Ed,o Pl°t i s l ar g e - An analysis of residuals shows that (Vd,m/K, m ) tends to decrease as the galaxy is 
redder, more luminous (massive), and of earlier type. The models allow us to explain the physics of these results, 
which imply a connexion between halo structure and luminous properties. The dynamical-to-baryon mass and 
dynamical mass-to-light (B and K) ratios at a given radius were also estimated. All these ratios, for observations 
and models, decrease with Sd,o (or SB) and do not correlate significantly with the galaxy scale, contrary to what 
has been reported in previous works, based on the analysis of rotation curve shapes. We discuss this difference 
and state the importance to solve the controversy on whether the dark and luminous contents in disk galaxies 
depend on SB or luminosity. The broad agreement between the models and observations presented here regarding 
the trends of the dynamical-to-baryon matter and mass-to-light ratios with several galaxy properties favors the 
ACDM scenario. However, the excess of dark matter inside the optical region of disk galaxies remains as the main 
difficulty. 

Key words. Cosmology: dark matter - Galaxies: disk - Galaxies: evolution - Galaxies: fundamental parameters - 
Galaxies: halos 



1. Introduction 

The information about the distribution of luminous and 
dark matter in disk galaxies, as well as on the correla- 
tions among the main parameters that characterize both 
components, offers a fundamental clue to understand how 
galaxies did form and evolve and what role did dark mat- 
ter play in these processes. Current models of disk galaxy 
formation and evolution, based on the hierarchical Cold 
Dark Matter (CDM) scenario, make certain predictions 
about these distributions and correlations (e.g., Firmani, 
Avila-Reese & Hernandez 1997; Dalcanton et al. 1997; Mo, 
Mao & White 1998; Avila-Reese, Firmani & Hernandez 
1998; van den Bosch 1998, 2000; Firmani & Avila-Reese 
2000; Avila-Reese & Firmani 2000). These predictions can 
be used to confront and interpret the available observa- 
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tions. This is the aim of the present paper, where we will 
compare some global properties and correlations of ob- 
served galaxies with the model ones. From the point of 
view of the models, the galaxy properties are the result of 
the combination of some fundamental (cosmological) pa- 
rameters, which follow continuous statistical distributions. 
Therefore, an observational sample as complete as possible 
in luminosities, surface brightnesses, morphological types 
and integral colors is crucial. 

More than 25 years ago, the interpretation of the ob- 
servations under the assumption of Newtonian dynamics, 
suggested the presence of dark matter in and around disk 
galaxies (Bosma 1978, 1981a,b; Rubin et al. 1980, 1982, 
1985; see also Rogstad & Shostak 1972). Two relevant 
(related) questions are whether the disk dominates or not 
the mass within some parts of the galaxies, and whether 
the dark halo is shallow or cuspy in the center (for re- 
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cent reviews see e.g., Bosnia 2002 and Salucci & Borricllo 
2001, respectively). The latter issue turned out crucial for 
testing the predictions of CDM models, and it is investi- 
gated currently by studying the inner rotation curves of 
observed dwarf and low surface brightness (LSB) galaxies, 
which are dark matter dominated systems. 

The rotation curve decomposition of individual galax- 
ies was commonly used to explore the fraction of dark 
matter within them (e.g., Carignan & Freeman 1985; van 
Albada & Sancisi 1986; Kent 1986). Besides of postulat- 
ing a stellar mass-to- light ratio at the observed band, 7a, 
this method requires to assume a dark halo mass distri- 
bution for the given galaxy. A popular technique is to fix 
7a (assumed constant with radius) at the highest possible 
value such that the "luminous" rotation curve component 
does not generate a higher rotation velocity than actu- 
ally observed in the inner part, the so called "maximum 
disk" hypothesis; the complement to the observed rota- 
tion curve is fitted to the assumed non-hollow dark halo 
model (typically a pseudo- isothermal sphere). For a maxi- 
mum disk model, the disk velocity component at its maxi- 
mum dominates over the halo. Luminous and high surface 
brightness (HSB) galaxies are typically well described by 
the maximum disk model, while this model for low lumi- 
nosity and LSB galaxies implies unrealistically high stellar 
mass-to-light ratios (see references in §3.4). The fractions 
of luminous and dark matter inferred with the individ- 
ual velocity decomposition methods are model-dependent 
and in certain circumstances the model may converge to 
an erroneous solution or to a non-reasonable values for 7a 
(Persic & Salucci 1991). Besides, for these methods de- 
tailed photometric and kinematic data are required. 

The fractional amount of dark matter in disk galax- 
ies have been also inferred by studying the kinematical 
properties of a set of rotation curves (Persic & Salucci 
1998, 1990a; Salucci & Persic 1999). The virtue of this 
method is that no assumption about a particular dark 
matter distribution should be done. Solely from the shape 
of the rotation curves and the definition of some radial 
scale connected to the luminous distribution (assumed ex- 
ponential) , it is possible to infer on average the fraction of 
dark and luminous matter at this scale. Persic & Salucci 
have found a dependence of the rotation curve slope mea- 
sured at 3.2 hd (/id is the total disk scale radius assumed 
equal to the optical scale radius) with the optical lumi- 
nosity of the galaxy, which translates into a dependence 
of the disk-to-total mass ratio at 3.2 on luminosity: 
Mdisk/Mtot oc Lg 4 . The samples used by these authors 
included only HSB highly-inclined late-type spirals with 
available Ha rotation curves and optical photometry. 

A major goal of the mass modeling method based on 
the shape of rotation curves was the inference of the so 
called "universal rotation curve" (Persic, Salucci & Stel 
1996). From about one thousand observed Ha rotation 
curves of late- type HSB spirals, synthetic rotation curves 
binned by luminosity intervals were generated. Assuming 
an exponential disk and a halo with a velocity profile pro- 
portional to r/(r 2 + a 2 ) 1 / 2 , where a is a core radius, the 



synthetic rotation curves could be well described by an 
universal profile, resulting from the quadratic sum of the 
disk and halo velocity components. The two parameters 
of this profile are the disk-to-total velocity ratio at 3.2fod 
and the core radius. Persic et al. found that, in order that 
the synthetic rotation curves could be fitted with the uni- 
versal profile, these parameters, in particular the former, 
should be a function of luminosity. Therefore, the univer- 
sal rotation curve is a function of luminosity in the sense 
that as the galaxy is less luminous, the larger is the slope 
of the rotation curve at 3.2h d . A direct implication of this 
result (which was in fact the inspiration for it) is that the 
dark-to- luminous mass ratio (at 3.2 /id) scales inversely 
with luminosity. 

In spite of the virtues of the rotation curve shape 
method to infer dark and luminous mass fractions, it is 
not free of uncertain assumptions (to be discussed in §4.2), 
which could bias the results. Moreover, other workers have 
shown several examples of rotation curves which deviate 
from the universal curve (Verheijen 1997; Bosma 1998), 
suggesting that besides the luminosity, other observational 
parameters, for example the surface brightness (SB), could 
be important, and therefore, the disk-to-total mass ratio 
could also depend on them. Thus, it is important to ex- 
plore alternative methods based on more complete sam- 
ples, in the sense of SBs, luminosities, and morphological 
types, to infer the dark and luminous fractions in disk 
galaxies. 

Some of these alternative methods are based on ob- 
served global galaxy parameters and relationships among 
them rather than in the local kinematics (e.g., Salucci, 
Ashman & Persic 1991). In the present paper, we follow 
this kind of methods, inferring from the observations (i) 
global parameters related to the amounts of luminous and 
dark matter in disk galaxies, and (ii) exploring how these 
parameters vary with galaxy properties. Part of the anal- 
ysis presented here is closely related to the fundamental 
plane of disk galaxies. However, a direct study of the fun- 
damental plane will be presented elsewhere. The results 
obtained will be compared with predictions of galaxy evo- 
lution models in the hierarchical clustering scenario, in or- 
der to test whether disks formed inside Cold Dark Matter 
(CDM) halos are realistic or not. 

A prediction of this scenario is that the luminous- 
to-dark mass ratio within the optical parts of the disks 
depends mainly on the disk surface density and on the 
global disk(baryon)-to-halo mass fraction, /d 1 (Firmani 
& Avila-Reese 2000; see also Dalcanton et al. 1997; Mo et 
al. 1998). A dependence with the mass (or luminosity) is 
not expected, unless astrophysical process, like feedback, 
introduce a strong dependence of /d on the mass. Recent 
calculations by van den Bosch (2002) show that the feed- 
back reduces /d in low massive galaxies but there is not a 
significant trend with the mass. In fact, in massive galax- 

1 We define / d = M haT /M v , where M bar is the total baryonic 
galaxy mass, referred some times also as the disk mass, Ma, 
and M v is the virial (total) halo mass. 
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ies not all the baryon fraction ends in the disk due to a 
large gas cooling time in these systems, so that both low 
mass and massive systems in average incorporate 0.4-0.6 
of the halo baryon fraction to the disk (van den Bosch 
2002), without any significant dependence on mass (i.e., 
/d is not expected to depend on mass). 

In spite of the great observational effort done in the last 
decades, there is actually a few number of homogeneous 
enough observational works reporting both photometric 
(optical and near-infrared) and kinematic information, in- 
cluding HI detection, for disk galaxies. In this paper, we 
present an extensive compilation from the literature for 
HSB and LSB disk galaxies with photometric parameters 
in both the B and K bands, rotation curves or veloc- 
ity line-widths and HI integral fluxes (see also Graham 
2002). After applying an uniforming procedure, the com- 
piled data is used to estimate several stellar and baryonic 
parameters of the disks. The most reliable data for this 
aim are in the near- infrared bands (H or K). The sample 
and the results in this and in a forthcoming paper can be 
used as a local reference for studies of the luminous and 
dynamical properties of disks in other environments and 
at higher redshifts. 

In §2, we present the compiled sample and describe 
the procedures to uniformizc the data and to estimate 
stellar and baryonic parameters from the observations. 
Sect. 3 is devoted to the study of the ratio of disk-to- 
total maximum velocity at their corresponding maxima, 
(^d,m/Vt, m )- This quantity is not defined at a given ra- 
dius and gives an approximate estimate of the dark and 
luminous fractions. For the kind of kinematical data we 
use here (equivalent line- width instead of detailed rotation 
curve), (Vd,m/Vt, m ) is the most direct quantity to com- 
pare with models, without introducing extra assumptions. 
In §§3. 1 simple composite halo/disk models predictions 
are presented, while in §§3.2 the description and results 
for semi-numerical galaxy evolution models in the ACDM 
scenario are given. In §§3.3 and 3.4, the (Vd, m /Vt, m ) ra- 
tio of observed and model galaxies and the dependence of 
this ratio on several galaxy parameters are presented and 
compared. In the light of the results obtained, the issue 
of maximum or sub-maximum disk is discussed. In §4, the 
dynamical-to-baryonic mass ratio (here instead of "total" 
we use the term "dynamical"), Md yn /Mb ar , and the mass- 
to-light ratios (both measured at some inner radius) for 
models and observations, as well as their dependences on 
several galaxy parameters are presented. For the obser- 
vational data, some assumptions should be done in order 
to estimate these ratios; we explore how much could affect 
these assumptions on our conclusions. In §4.2, we compare 
our results with those in previous works and discuss the 
differences. Finally, a summarizing discussion is given in 
§5. 

2. The Sample 

In order to trace the baryonic and dark matter content 
of disk galaxies it is important to explore a wide range of 



galaxy types, magnitudes, sizes, and surface brightnesses 
as possible. Furthermore, high-quality surface photometry 
in the optical and near-IR pass-band as well as informa- 
tion about the rotation curve or at least the HI or Ha line- 
width, and the total HI gas flux are needed. After an ex- 
tensive search in the literature, we surprisingly have found 
only a few (homogeneous enough) observational sources to 
compile a sample of disk galaxies with the required data 
(see also Graham 2002). From these sources, three main 
sub-samples of HSB and LSB galaxies were collected: 

(1) de Jong sub-sample (de Jong & van der Kruit 1994, 
and de Jong 1996a) reports bidimensional B, V, R, I and 
A-band photometry for 86 undisturbed field spirals se- 
lected from the Uppsala General catalogue (UGC). Only 
galaxies classified equal or later than SI or SB1 and with 
an inclination limit < 51° were included. This is a diam- 
eter limited sample (UGC is expected to be complete up 
to 2 arcmin red diameters) that can be transformed into 
a volume limited sample. 

(2) Verheijen (1997) sub-sample (see also Verheijen & 
Sancisi 2001) is a selection of galaxies from the Ursa Major 
Cluster. This is one of the least massive and most spiral- 
rich nearby cluster. The sample is complete for galaxies 
later than Sab to a limiting apparent Zwicky magnitude 
mzw = 15.2, although it includes galaxies fainter than 
that cutt-off. It contains B,V,R, J, A'-band photometry 
and kinematic (HI observations) for HSB and LSB disk 
galaxies. The kinematic data were obtained by the authors 
only for galaxies (52) with inclinations greater than 45° 
and include synthetic rotation curves as well as equivalent 
line- widths. Notice that although it is suspected that clus- 
ter environment affect the properties of the constituent 
galaxies, in the case of the Ursa Major, the cluster may 
be so young that it's members are more representative of 
a field population. 

(3) Bell et al. (2000) sub-sample, was obtained par- 
tially by themselves and compiled partially from the lit- 
erature, mainly from de Blok, van der Hulst & Bothun 
(1995) and de Blok, McGaugh & van der Hulst (1996), 
with the condition to have blue central surface brightness 
Ms,o > 22. 5magarcsec~ 2 and diameters at the 25 mag 
arcsec~ 2 isophote larger than 16". Their sample consist 
of 26 LSB galaxies with available bidimensional B— and 
AT— band photometry and an inclination limit < 61°. The 
sample is by no means complete but it is designed instead 
to span as wide a range as possible of observed LSB galaxy 
parameters. Notice that although LSB galaxies follow the 
spatial distribution of HSB galaxies, they tend to be more 
isolated from their nearest neighbors than HSB galaxies 
(c.f. Bothun et al. 1993). 

The Milky Way and Andromeda galaxies were also 
included so that a direct comparison of their properties 
w.r.t. the galaxies in our sample could be inferred. 

The galaxies from the Verheijen sub-sample and some 
LSB galaxies from the Bell et al. sub-sample have mea- 
sured synthetic HI rotation curves. Thus, their maximum 
total velocity rotation, Vt.m, is known. For the rest of the 
galaxies in our collected sample, only equivalent velocity 
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line-widths are reported. In order to uniformize the data 
as much as possible, we used the velocity line-width W20 
(defined at the 20% level) for all the galaxies in the sam- 
ple. Only a few (LSB) galaxies in the sample lack of W20 
but have measured rotation curves; in these cases we have 
used the Vt, m inferred directly from the rotation curve. If 
properly corrected, W20 is « 2Vt, m as Verheijen (1997, p. 
197) have shown (see also Verheijen & Sancisi 2001). Thus, 
the dynamical information available for our sample will be 
Vt, m - The radius where V t peaks, symbolized here as r poa k, 
is not obviously related to a typical radius of the disk but 
it can be estimated as a multiple of hd (see §3). For the 
Verheijen sub-sample, we adopt his W20 values, while for 
the de Jong sub-sample and most of the LSB galaxies from 
Bell et al. sub-sample, W20 was extracted from the Lyon 
Extragalactic Database (LED A) database 2 . The integral 
fluxes in the 21-cm line (used to calculate the neutral gas 
masses in galaxies) were taken from Verheijen & Sancisi 
(2001), LEDA, and dc Blok et al. (1996). 

A condition imposed to our sample is that galaxies 
should be in a restricted range of inclinations (35° < 
i < 80°). Otherwise either the surface brightness (here- 
after SB) profiles or the rotation velocities are less reli- 
able. The inclinations were calculated from the photomet- 
ric minor-to-major semi-axis ratios (typically in a near- 
infrared band) reported in the source papers and assum- 
ing a fixed thickness parameter q = 0.2. We exclude from 
the sample galaxies with clear signs of interaction (mainly 
from the Verjeihen sub-sample) and with rotation curves 
which are still increasing at the last measured outer point. 
Only some LSB galaxies showed this feature. Fortunately, 
there are published synthetic rotation curves for most of 
the LSB galaxies in our sample. The final sample consists 
of 78 galaxies: 42 out of 86, 29 out of 52, and 5 out of 23 
from the de Jong, Verheijen and Bell et al. sub-samples, 
respectively, plus the MW and Andromeda galaxies. The 
Bell et al. sample contains actually 26 LSB galaxies, but 
three of them were taken from the de Jong sample. 

The local distances to the galaxies were calculated us- 
ing the kinematic distance modulus given in LEDA. The 
value of the Hubble constant used in this database is 
the same we assume here, Ho = 70kms~ 1 Mpc~ 1 . Since 
the LSB galaxies from the Bell et al. sub-sample arc 
not included in LEDA, their distances were taken di- 
rectly from Bell et al. (2000) correcting them for the 
H value used here. From the given distances, the corre- 
sponding rcdshifts were obtained. Using these redshifts, 
we then calculated the luminous distances as: Dl = 

[l + i(l — qo)z] , where q is the deceleration parame- 
ter and c is the speed of light. For the cosmology used in 
this paper (Q m — 0.3,17a = 0.7), qo = —0.55. The above 
equation is a good approximation for small redshifts. 

Merits and limitations of the sample.- The sample is 
not complete in any sense, but it is representative of the 
basic parameters required: morphological type, luminos- 
ity, SB and color. In Fig. 1, the number distribution of 

2 http://leda.univ-lyonl.fr 
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Fig. 1. Number distributions of morphological types, 
B— band corrected absolute magnitudes and SBs, and cor- 
rected (B — K) colors for the recompiled sample of 78 nor- 
mal HSB and LSB disk galaxies. The contribution of the 
HSB galaxies is shaded. 



morphological types, corrected B— band magnitudes and 
SBs, and corrected (B — K) colors are presented (for the 
corrections, see §2.1). The number distribution of types 
is similar to that in other samples of nearby galaxies 
(c.f., Hernandez-Toledo et al. 1999; Jansen et al. 2000), 
with a high peak at T = 4, 5. The distribution of abso- 
lute magnitudes indicates that the sample is representa- 
tive of the local galaxy luminosity function in the bright 
part (M a b s ^ — 19.5) but lacks enough representation for 
fainter luminosities. In fact, the faint part of the luminos- 
ity function consists mostly of dwarf galaxies. The physics 
of this group of galaxies is different from that of normal 
disk galaxies. The aim of the present paper is to investi- 
gate only the latter ones. The SB distribution presents a 
bimodality (Tully & Verheijen 1997). The sample is not 
complete at very LSBs. This is consistent with the lack of 
faint galaxies in the magnitude number distribution: very 
LSB galaxies are typically dwarf galaxies. 

Since we adopt the raw disk central SB and scale ra- 
dius (^a,o and h\) directly from the above source samples, 
it is important to check that the the various methods used 
to infer these parameters do not introduce systematic dif- 
ferences in our compiled sample data. For the de Jong 
sample, a two-dimensional exponential bulge-bar-disc de- 
composition was used (de Jong 1996). Bell ct al (2000) 
used a one-dimensional bulge-disc decomposition (with ei- 
ther an exponential or an r 1 / 4 law bulge). Verheijen used a 
'marking the disc' fit, where the contribution of the bulge 
to the disc parameters is minimized by visually assess- 
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ing where the contribution from the bulge component is 
negligible. Dc Jong (1996a) showed that the one- and two- 
dimensional methods give comparable results to typically 
better than 20% in terms of surface brightness (SB here- 
after) and 10% in terms of scalelengths, with little signif- 
icant bias. 

After a comparison of the de Jong and Verheijen sam- 
ples in the L\-h\ diagram for the B and if bands, we 
found no systematic differences. Notice that the selection 
criteria applied in de Jong(1996a), implies that small LSB 
galaxies are excluded. This selection effect could intro- 
duce a systematic difference with the Verjeihen galaxies. 
However, only a few galaxies from the later author fall 
below the selection limit of de Jong in the ^Kfl-hx dia- 
gram (Ho = 70kms~ 1 Mpc~ 1 ; this value is used throughout 
the paper). We also checked that the scale radii reported 
by these authors do not correlate with inclination (e.g., 
Giovanelli et al. 1984). When the bulge is not taken into 
account in the fitting, such a correlation could appear, but 
that is not the case for the Verheijen data. 

In conclusion, we do not find significant systematic dif- 
ferences in /za,o and h\ among the galaxy sub-samples due 
to the different profile fitting procedures used. This is con- 
sistent with de Jong (1996a) conclusion mentioned above. 
In fact, the dominant source of error in the determination 
of the galaxy photometric parameters is due to the uncer- 
tainty in the sky level (Bell et al 2000; de Jong 1996a). 

2.1. Data Corrections 

The total magnitudes were corrected for Galactic extinc- 
tion (Schlegel, Finkbeiner & Davis 1998), if— correction 
(Poggianti 1997), and internal extinction. Since the frac- 
tion of dust seems to be larger for bigger galaxies 
according to empirical (e.g., Wang & Heckman 1996; 
Tully et al. 1998) and theoretical (e.g., Shustov, Wicbe 
& Tutukov 1997) arguments, we used the empirical 
velocity- (luminosity-) dependent extinction coefficients 
determined by Tully et al. (1998) for a large sample of 
galaxies (HSB and LSB) in clusters. The if'-band data 
from the Verheijen sub-sample were transformed into K- 
band adopting if' — if = 0.19(if — if) (Wainscoat & Cowie 
1992) and the average (H—K) color of each morphological 
type given in de Jong (1996c). 

The SBs were corrected for Galactic extinction, 
if— correction, cosmological SB dimming, and inclination 
(geometrical and extinction effects). For the latter cor- 
rection, we followed Verheijen (1997), considering LSB 
galaxies as optically thin in all bands. The HSB galax- 
ies were also considered as optically thin in the if-band. 
We have defined LSB galaxies as those whose disk central 
SB in the if— band after SB correction is larger than 18.5 
mag/arcsec 2 (Verheijen 1997). 

The 21 cm line-widths at the 20% level, W 20 , for the 
Verheijen susbsample were taken directly from Verheijen 
& Sancisi (2001). Galaxies with W20 taken from LEDA, 
were de-corrected (for instrumental corrections) to get the 



raw data (Paturel et al. 1997) and then, corrected again for 
broadening due to turbulent motions and for inclination, 
this time following Verheijen & Sancisi (2001). 

The final sample data are presented in Table 1. The 
name and morphological type are given in Col. (1). The lu- 
minous distance and the galaxy inclination are presented 
in Col. (2). Raw apparent magnitudes and corrected lumi- 
nosities in the B and if bands appear in Cols. (3) and (4), 
respectively. The solar magnitudes used were Mbq = 5.48 
and Mkq = 3.41. The raw and corrected central SBs 
in the B— and if— bands appear in Cols. (5) and (6), 
while Col. (8) shows the corresponding disk scale radii 
(in kpc). In Col. (7), the corrected integral and central 
(B — if) colors are given. Finally, Col. (9) gives the cor- 
rected line-widths at the 20% level, W20, and the HI mass, 
M H i(M Q ) = 2.36 x 10 5 D|(Mpc) j* ^(Jykms" 1 )^. 

2.2. Composite quantities 

The stellar mass and disk central surface density, M s and 
S S) o, are derived from the if —band luminosity (which in- 
cludes the bulge) and central SB, respectively. It seems 
that, even in the if band, the stellar mass-to-light ratio, 
7k, is not constant. The application of population syn- 
thesis techniques to simplified galaxy evolution models for 
HSB galaxies show that jk depends mainly on the inte- 
gral color (Tinsley 1981; Bruzual 1983; Bell & de Jong 
2001). From Table 1 (model "formation epoch with bursts 
(vi)") in Bell & de Jong (2001), we infer 



Log 7K = -0.91 + 0.21(B - if). 



(1) 



For low SBs this relation breaks down. The various ro- 
tation curve decomposition models applied by Verheijen 
(1997) to his sample of LSB galaxies, actually shows that 
jk increases toward bluer galaxies. A linear eye fit to the 
jk — (B — if) relation of his LSB galaxies, assuming a 
Hcrnquist halo model and using his constrained decompo- 
sition method, gives 



j K = 2.75 - 0.75(5 -if'), (B -if') <3 
1K = 0.5, (B - if') > 3. 



(2) 



This result is similar to that inferred from other decom- 
position methods and for the pseudo-isothermal halo used 
by Verheijen (1997; see Fig. 10 in Chapter 6). The scale 
radius of the stellar disk, h s , is assumed here to be equal 
to the scale radius in the band if, Hk- Most of the stars 
in the disk are traced in this band. 

The disk gas mass, M g , is estimated as 



M„ = 1.4M H i 



1 + 



M H2 



(3) 



where the factor 1.4 takes into account helium, and Mh 2 
is the molecular hydrogen mass. The M# 2 /M#j ratio has 
been found to depend on the morphological type T (Young 
& Knesek 1989). From the latter paper, McGaugh & de 
Blok (1997) estimated that M H2 /M H i = 3.7 - 0.8T + 
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0.043T 2 . For T< 2 this empirical fitting formula overesti- 
mates strongly the gas mass in galaxies, thus we assume 
M H2 = 0. 

The estimation of the gas and baryonic galaxy param- 
eters is important to understand several structural and 
dynamical aspects of galaxies and to compare with the 
theoretical predictions. The galaxy baryonic mass is Afbar 
= M s + M g , and the gas mass fraction is defined as / g 
= Mg/Mbar- Unfortunately, we do not have information 
about the gas surface density profiles (S a ,o and h g ) for all 
galaxies in our sample. We need them in order to calcu- 
late the baryonic disk central surface densities and scale 
radii. The observed density profiles of neutral gas in disk 
galaxies are diverse and there seems to exist a character- 
istic profile for each morphological type (Broeils & van 
Woerden 1994). The HI distribution in the center of Sb- 
Scd galaxies is typically flat or even has a hollow, but 
the H2 distribution for these galaxies rapidly increases 
to the center compensating the HI deficiency. For early- 
type galaxies, the HI profile is roughly exponential and 
H2 is deficient. In a very rough approximation, the to- 
tal gas distribution (HI+H2) in disk galaxies can be de- 
scribed by an exponential profile with scale radii 2.0-3.0 
times larger than the stellar scale radii in the near-infrared 
bands (e.g., Olivier, Blumenthal & Primack 1991; Marsh 
& Helou 1995). We will assume that the total gas surface 
density follows an exponential distribution with an scale 
radius 3 times Kk ■ Thus, £ Sj o = M g /2ir(3hK) 2 ■ The bary- 
onic disk central surface density is then calculated as £d,o 
= T, S fi + S s ,o- The corresponding baryonic disk scale ra- 
dius will be h d = /ix[(S Sj o + 9£ gi o)/5]d,o)] ' 5 - bi fact, our 
main results in §3 will not be affected by including the gas 
disk parameters, but some influence is expected. This is 
why we tried to use the information we have at least about 
the global gas content, by introducing and "effective" gas 
scale radius and surface density. 

In Table 2 we present the composite quantities de- 
scribed above for the sample galaxies in Table 1. Cols. (2) 
and (3) contain the disk stellar mass and central surface 
density (M s and £ s ,o)- Cols. (4) and (5) the total disk gas 
mass and mass gas fraction (M g and / g ). Finally, Cols. (6), 
(7) and (8) list disk baryonic mass, central surface density 
and scale radius (Mb ar , ^d.cij and ^d, respectively). 

3. Disk-to-total maximum velocity ratio 

We are interested in exploring the mass amounts of lu- 
minous and dark matter in disk galaxies and how these 
amounts do depend on different galaxy properties. Our 
goal is to compare the observational results with theo- 
retical predictions. As mentioned in the introduction, the 
analysis of the rotation curve shapes offers a powerful 
method to estimate the fraction of luminous and dark 
matter in galaxies at a given radius (3.2/id, for example). 
The approach used here is conceptually different. In our 
sample, we have information only about the HI velocity 
line- width, W20, which approximately corresponds to twice 
the maximum total rotation velocity Vt,m (see introduc- 



tion), but not about the rotation curve shapes. On the 
other hand, with the information about the disk struc- 
ture that we have (assuming exponential stellar and gas 
distributions), the velocity component due to the to- 
tal disk can be calculated (eq. El below) , for example at 
the radius where V ( i attains its maximum, 2.2h d . Thus, 
without introducing further assumptions, we may define 
the ratio of disk maximum velocity to total (or dynami- 
cal) maximum velocity, (Vd,m/Vt. m ), which can be directly 
compared with theoretical predictions. 

The (Vd.m/Vt^m) ratio is not defined at a given radius 
but it can be related to the dynamical-to-baryon (disk) 
mass ratio, Md yn /Afbar, defined at the radius where the 
total rotation curve peaks, rpeak. Assuming an exponential 
disk (Mbar oc £d,o^d 2 ), one obtains: 

/ M dyn \ ^ Vt.nApcak/G' ^ Vt,m 2 xfod 

VM bar / P eak f L 2nY, dfi h d 2 fL(x)T; dfi h d 2 

where r poa k = xh d , fa is the fraction of the total disk mass 
at r poa k, and Vd, m is calculated according to eq. © below. 
So, in order to calculate M dyn /Mb al , one needs to define 
the coefficients x and /l(x); they actually may change 
from galaxy to galaxy, depending on the disk SB, lumi- 
nosity, and halo density profile (x « 2.2 and Jl ~ 0.64 
for disk dominated galaxies, and probably, x > 2.2 and 
/l > 0.64 for halo dominated galaxies). Previous works, 
based on the analysis of rotation curves, have shown that 
Mdy n /Mbar (at 3.2h d ) has a strong dependence on lumi- 
nosity (Persic & Salucci 1988,1990a; Persic et al. 1996). In 
§4 we will present results for Md yn /Afb a r- Here we concen- 
trate on (Vd, m /Vt tm ), which is more directly comparable 
to model results than M dyn /Mi, ar , and may offer a first 
idea of how the fractions of dark and luminous matter in 
galaxies do depend on several galaxy properties. 

In order to calculate Vd,m from the observations, we 
use the expression of for a flattened exponential disk 
(Freeman 1970): 

V*(y)=4irGZ lifi h d y 2 I(y), (5) 

where y — r/(2h d ), G is the gravitational constant, and 
I(y) is a function composed of modified Bessel functions 
of the first and second kinds. The function (JSJ has its max- 
imum at 2.2h d . Notice that £d,o and h d are not directly 
observed quantities (see §2.2). Although the assumptions 
made to calculate them from the observational data are 
robust, we will check a posteriori whether our results are 
significantly dependent on these assumptions. 

According to Sacket (1997, and more references 
therein), the maximum disk hypothesis applied to individ- 
ual Sb-Sc galaxies implies a baryonic disk mass (includ- 
ing gas) which provides (85 ± 10)% of the total rotation 
support at 2.2h d (2.2\ik)- The lower limit is occupied by 
galaxies with large bulge; only in these cases the bulge 
contribution to Vt im may become significant. From the 
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analysis of rotation curve shapes and assuming the max- 
imum disk hypothesis, Persic & Salucci (1990b) reported 
-Mbar/Mdyn ratios at 3.2/id increasing with L# from 0.3 
to 0.8. This implies that (Vd/Vt) 3 2 ^ d oscillates between 
0.57 and 0.91. The Vd/Vt ratio increases from 3.2 to 2.2hd; 
by using the rotation curve decompositions in Fig. 3 of 
Persic et al. (1996), we estimate that (Vd/V t ) 2 2 f ld would 
oscillate between 0.75 and 0.95, in excellent agreement 
with Sacket (1997). Although (Vd <m /V t , m ) is not formally 
equal to (Vd/V t ) 2 2 ^ d , as a reference we will say that the 
maximum hypothesis implies (Vd, m /Vt :m ) » 0.85 ± 0.1. 
Since the individual rotation curve decomposition meth- 
ods usually take into account the disk thickness (q > 0), 
we will correct Vd, m as given in eq. © for disk thickness. 
Assuming a minor-to-major disk axis ratio q of 0.2, Vd, m 
becomes approximately 5% lower than for the infinitely 
thin disk case (e.g., Burlak, Gubina & Tyurina 1997). We 
will apply this correction factor to the disk+halo compos- 
ite mass models presented below as well to the observa- 
tional data. 

Before studying and comparing the inferred 
(Vd,m/Vt, m ) ratios with results from galaxy evolu- 
tion models, we shall explore how this ratio is expected 
to depend on different disk and halo parameters. To this 
purpose, simple disk+halo composite mass models will 
be used, and then they will be compared to the complex 
evolutionary models. 

3.1. Composite mass models using analytical halo 
profiles 

Here we will analyze the rotation curve decomposition cor- 
responding to exponential disks in centrifugal equilibrium 
inside pseudo- isothermal (e.g., de Blok & McGaugh 1997) 
and Navarro, Frenk & White (1997; hereafter NFW) dark 
halos. A description of how the halos and disks are con- 
structed is given in Appendix A. In Fig. 2a we present 
(Vd, m /Vt,m) vs. the central disk surface density, £d,o> f° r 
the virial (halo) mass M v = 3.5 10 11 M Q and for the two 
halo models. Upper and lower curves correspond to the 
pseudo-isothermal and cosmological NFW halos, respec- 
tively. In the former case, the assumed core radius was 3.8 
kpc (see Appendix A). The solid lines are for the global 
baryon-to-halo mass ratio /d= 0.05, while dashed lines 
are for / d = 0.08 (upper) and 0.03 (lower). For the NFW 
halos, we have used the average concentration parameter 
for a given mass taken from Eke, Navarro & Steinmetz 
(2001). The ACDM cosmology with fl m = 0.3, Q A = 0.7, 
h = 0.7 is used here. The disk central surface density is 
determined mainly by the halo spin parameter A, but de- 
pends also on M v and fd (see Appendix A). Figure 2 is 
for a range of A from 0.015 to 0.2. The trend of decreas- 
ing (Vd, m /Vt ;m ) as £d,o decreases is clearly seen for both 
halo models. Obviously, for the pseudo-isothermal halos, 
the disk velocity component dominates more than for the 
cuspy NFW halos. 
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Fig. 2. (a). Correlation between (Vd, m /Vt, m ) and Sd,o for 
disks in centrifugal equilibrium inside a pseudo-isothermal 
(upper curves) and a NFW (lower curves) halo of M v = 
3.5 10 11 M Q . See details about the model parameters in the 
text. Solid lines are for / d = 0.05, while upper and lower 
dashed lines around the solid curves are for /d=0.08 and 
0.03, respectively. The approximate region of maximum 
disk (see text) is bounded by the dot-dashed lines, (b). 
Same as in panel (a), but for the NFW model, taking into 
account the disk gravitational dragging over the halo. The 
shaded region encompasses de ±lcr scatter in concentra- 
tions around the central model. The solid upper and lower 
curves are for M v = 3.5 10 12 M Q and 3.5 1O 1O M with the 
corresponding average concentration parameters and /d= 
0.05. (c) and (d). Relation between (Vd, m /Vt,m) and the 
disk(baryon)-to-dynamical mass ratios inside 2.2/id and 
5/id for the same models shown in the left panels. In panel 
(c) thick and thin lines are for the NFW and pseudo- 
isothermal halos, respectively. Botton in panel (d), the 
linear fittings (in Log-Log) arc given. 

While the pseudo-isothermal halo profile is an adhoc 
model introduced to fit some observed rotation curves, the 
NFW profile is a theoretical prediction for the structure 
of the halos in the CDM cosmogony. In the latter case, it 
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is important to take into account that the inner structure 
of these halos is affected by the formation of the bary- 
onic galaxy, which gravitationally contracts the original 
halo. The usual way to calculate this effect is by assuming 
a radial adiabatic invariance during the formation of the 
disk (e.g., Flores et al. 1993; Mo et al. 1998; see Appendix 
A). In Fig. 2b, the same curves for the NFW model pre- 
sented in panel (a) are shown, but this time taking into 
account the halo contraction due to disk formation. Two 
more masses, M v = 3.5 I0 lo M Q and 3.5 1O 12 M are added 
in this panel. As one sees, when the halo model is origi- 
nally of NFW type, after disk formation, even disks with 
high surface densities tend to be of sub-maximal type: 
(Vd, m /Vt,m) almost never exceeds ~ 0.75 — 0.80 3 . The ef- 
fects of the halo contraction on the (Vd, m /Vt im ) ratio are 
non-significant only for the lowest surface density disks. 
For the M v = 3.5 lO n Af model with / d =0.05, we also 
show the scatter produced in the (Vd )m /Vt im )-£d.o dia- 
gram by the different halo concentrations. The shaded re- 
gion encompasses ±ler deviations in c (AlogCNFW ~ 0.18, 
Bullock et al. 2001; the lower limit corresponding to the 
+ 1(7 value). 

A comparison in the (Vd )m /Vt,m)-£d,o P^ ane °^ ^ ne 
pseudo-isothermal and NFW halos (taking into account 
the disk gravitational drag in the last case), shows that 
as Sd,o is higher, the (Vd, m /Vt, m ) ratio for the pseudo- 
isothermal models becomes larger than that for the con- 
tracted NFW models. Thus, it is at high surface densities 
(SBs) where the the differences in the (Vd, m /Vt,m) ratio 
are more sensitive to whether the halo is cuspy or shal- 
low. On the other hand, neither the concentration nor the 
mass introduce significant differences in (Vd im /14 im )- 

For the NFW models with M v = 3.5 1O 1O M (upper 
thin solid line) and 3.5 10 12 M Q (lower thin solid line), 
we also have used the average NFW concentrations corre- 
sponding to these masses (Eke et al. 2001) and /d = 0.05. 
Massive systems result with lower (Vd, m /Vt, m ) ratios than 
the less massive ones, which apparently would imply that 
the former are more dark matter dominated than the lat- 
ter. In fact, the average concentration for less massive ha- 
los is higher than for the more massive ones. However, it 
is important to have in mind that the parameters of the 
modeled disks are connected to the scale of the halo: h^ 
and Ed,o depend on M v (see Appendix A). For a given A 
(constant Sd,o), oc , with [3 « 1.5, i.e. shallower 

than 2 (Firmani & Avila-Reese 2000, hereafter FA00). 
For constant Xd,o> (Vd, m /Vt, m ) oc /id°' 5 /Vt, m . Therefore, 
(Vd,m/V t , m ) oc V t , m ~ ' 25 , i.e. (Vd,m/14,m) decreases with 
increasing 14, m ( or M v ). This effect is partially compen- 
sated by the fact that more massive halos are less con- 
centrated (lower (Vd, m /Vt jm )). The final result is a small 
dependence of (Vd. m /Vt im ) on mass. For a variation of 2 

3 As mentioned above, the Vd. m calculated using eq. JSJ is 
lowered by 5% in order to take into account for the disk thick- 
ness (q = 0.2). This correction enters also in the total rotation 
velocity because this last is the sum in quadrature of the disk 
and halo velocity components. 



orders of magnitude in mass, the change in (Vd, m /Vt,m) 
is around 16% for a medium value of £d,o = 400M Q /pc 2 
(Fig. 2b). The change of (Vd im /14,m) f° r a ±lf in concen- 
tration is close to 12% for the same central surface density. 

The dependence of the (Vd, m /14,m) ratio on /d is also 
small. Curves with dashed lines in Fig. 2b are for the same 
central model (M v = 3.5 10 n M Q and average concentra- 
tion parameter) but /d from 0.08 (upper curve) to /d= 
0.03 (lower curve). The (Vd, m /Vt, m ) ratio does not change 
more than 10% within this range of /d for £d,o around 
400M Q /pc 2 . 

The (Vd, m /14,m) ratio is directly related to the disk 
and halo mass fractions. In order to gain a more quantita- 
tive feeling of the (V^m/^.m) ratio, Fig. 2, panels (c) and 
(d) show the corresponding disk-to-dynamical mass ratios 
at 2.2/id and 5/id for the same curves shown in panels (a) 
and (b). When these ratios are 0.5, the disk and the halo 
masses are equal at 2.2/id and 5/id, respectively. The rota- 
tion curve peaks typically at 2.2/id, while bhd corresponds 
approximately to the optical disk radius. For the gravita- 
tionaly contracted NFW halo model (panel d), the maxi- 
mum disk case ((Vd,m/Vt,m) > 0-?5) implies that 53% or 
more of the dynamical mass at 2.2/id should be in the disk. 
At 5/id this lower limit decreases down to 40%. Such 
high values can be attained only for ultra-high surface 
density disks. It is interesting to note that, for the NFW 
halos with contraction and for a large range of disk masses 
and surface densities, (Vd, m /Vt,m) ~ (M dyn /M har ) 2 ^ = 

(v<*/v t ) 2 . 2/ld . 

3.2. Self-consistent galaxy evolution models 

In previous papers (Avila-Reese et al. 1998; FA00; Avila- 
Reese & Firmani 2000; see also van den Bosch 1998, 2000, 
2002) semi-numerical models of disk galaxy formation and 
evolution in the hierarchical ACDM scenario were devel- 
oped. These models include self-consistently halo forma- 
tion and evolution, disk star formation (SF) and feedback 
processes, the gravitational dragging of the halo due to 
disk formation, secular bulge formation and other evolu- 
tion processes. The initial and border conditions of the 
models are M v , /d, the halo mass aggregation history 
(MAH) and the spin parameter A. A short description 
of the main ingredients of the models is presented in 
Appendix B. 

Before analyzing the predictions of these models re- 
garding the (Vd, m /Vt,m) ratio, it is important to remark 
the following three effects, (i) The halo contraction due 
to disk formation in the models was calculated for ellip- 
tical orbits, i.e. the adiabatic invariance was not limited 
to only radial orbits, as was done in the previous sub- 
section. Particles within the virialized ACDM halos have 
orbits with ratios of apapsis-to-periapsis around 0.15-0.40 
(Ghigna et al. 1998). Therefore, the adiabatic halo con- 
traction for the evolutionary models presented here is of 
lower amplitude than in the case of assuming radial or- 
bits (previous subsection), (ii) The disks formed within 
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ACDM halos, assuming detailed angular momentum con- 
servation, have actually some excess in the center (cusp) 
w.r.t. the exponential surface density distribution (FAOO; 
Avila-Reese & Firmani 2002; see also Bullock et al. 2001, 
van den Bosch 2002), which give rise to a bulge formed 
secularly (Avila-Reese & Firmani 2000). The effects (i) 
and (ii) imply that Vd, m and the (Vd, m /Vt im ) ratio mea- 
sured in our evolutionary models should be larger than 
in the simple composite mass models, where the orbits 
were assumed radial and the disk surface density expo- 
nential (see previous subsection), (iii) Our model disks 
are not infinitely thin; they have a vertical density and 
velocity structure determined by the SF, feedback and in- 
fall mechanisms included in the modelation. The gravi- 
tational potential of the disk is calculated taking into ac- 
count its three-dimensional structure. A comparison of the 
Vd,m measured in the model galaxies with the Vd. m calcu- 
lated according to eq. (|SJ) shows that the former is on 
average ~ 7% smaller than the latter. In fact we have cor- 
rected (uniformly) the composite mass models (and the 
observational data) to take into account the disk thick- 
ness. The 5% correction in Vd, m calculated according to 
Burlak et a. (1997) is close to the 7% measured on average 
in the evolutionary models. 

In Fig. 3a we present the results for 63 models grouped 
around three halo masses, M v ss 3.5 10 10 ,3.5 10 11 , and 
3.5 10 12 Mq (triangles, squares and circles, respectively) 
with MAHs and A's obtained randomly from their cor- 
responding (cosmological) statistical distributions, and 
/ d = 0.05 (FAOO; Appendix B). The solid curves are the 
same curves plotted in Fig. 3b for the three masses, the 
average concentration parameters, and f^— 0.05. The evo- 
lutionary models seem to lie slightly higher than the cor- 
responding composite mass models in the (V^m/K.m) vs. 
Sd.o plane, in particular at ultra-high surface densities, 
where the disk cusps are more pronounced. Although an 
attempt to introduce a correction due to disk thickness 
has been done for the composite mass models, it is im- 
portant to keep in mind the other two effects mentioned 
above as well several evolution processes included in the 
modelation. Nevertheless, the simple disk+halo composite 
mass models offer a reasonable description of results from 
the evolutionary models and deserve to be included. These 
models can be easily reproduced by other workers in the 
field interested in using these kind of theoretical results. 

The model (Kl,m/^t,m) ratio is determined mainly by 
the spin parameter A (FAOO), therefore (Vd,m/Vt,m) cor- 
relates mainly with Ed,o as seen in Fig. 3a. The models 
show also a slight dependence of the (Vd, m /Vt, m ) ratio on 
mass (see the segregation of triangles, squares and circles 
in Fig. 3a). The reason of this dependence has been al- 
ready explained in the previous subsection. For a given 
mass, the scatter seen in Fig. 3a is due to the scatter in 
the halo structures. This scatter has its origin in the dis- 
persion of the MAHs (Avila-Reese et al. 1998; Wechsler 
et al. 2002). For perfect NFW halos, this scatter is com- 
pletely described by the scatter in the NFW concentra- 
tions. Thus, more concentrated halos and more massive 
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Fig. 3. (a). Same as in Fig. 2 but for observations 
(empty circles) and evolutionary model results (solid 
triangles, squares, and circles corresponding to M v = 
3.5 10 10 M Q ,3.5 10 n M Q , and 3.5 10 12 M Q , respectively). 
The models are for /d= 0.05. Four galaxies, which are out- 
liers in this plot, as well as the Milky Way and Andromeda 
galaxies are indicated by their names. For indication, the 
composite mass models presented in Fig. 2 for the pseudo- 
isothermal halo (upper curves) and the NFW halo with 
gravitational drag (lower curves) are also shown. The solid 
lines are for M v — 3.5 10 n M© while the dashed ones 
are for 3.5 10 10 M Q (upper curves) and 3.5 1O 12 M (lower 
curves), (b). Means per Sd,o bins and their correspond- 
ing standard deviations for the observational data (thick 
dashed line and shaded area) and the model results (thick 
solid line and empty frame) showed in panel (a). The ar- 
rows from the thick curve are an estimate of how much 
the models would decrease if fa is decreased to 0.03. 



galaxies tend to have lower (Vd, m /^t,m) ratios for a given 
Ed,o- The differences are actually small. 
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3.3. Observational data 

In Fig. 3a we plot the (Vd,m/H,m) ratio vs. Sd,o inferred 
from the compiled observations in §2. The scatter is large 
and before comparing to models, we shall explore whether 
the observed correlation is significant and direct. In Fig. 
3a, four galaxies have been identified as outliers (shown 
with their names in the plots). These galaxies were ex- 
cluded from any statistical analysis. 

In a first exploratory phase, we look for any correla- 
tion of the (Vd, m /Vt,m) ratio (Log(V d , m /Vt,m)) with sev- 
eral galaxy properties (independent variables): Log£d,0; 
Log/id (or LogMbar), (B — K) color, morphological type 
T, and Log/ g . We applied backward step-wise multiple 
linear regressions, where F-tests are used to compute the 
significance of each independent. The elimination of non- 
significant variables leads us to a bidimensional linear cor- 
relation, where the two more significant independents are 
Sd,o and T, although it could also be S dj o and (B — K). 
The correlation is: 



i\r l\r \ „ v 0.20±0.02r r 0. OliO. 003 



(6) 



with the intercept equal to io~°- 72±0 - 06 and an ad- 
justed correlation coefficient R=0.79 (Sd,o is in unities 
of M Q pc~ 2 ). The dominant correlation is with £d,o- 
Notice that although a direct linear regression ren- 
ders (Vd,m/% m ) tx £ d ,o ' 14±0 ' 02 , with an intercept of 
1(r o.54±o.04 and r-o.75, as is seen in eq. ©, the cor- 
relation significance is improved slightly introducing T, 
although the slope for T is negligible. The reason for this 
is that E d ,o and T are not completely independent among 
them; £d,o (anti)correlates slightly with T. 

The most important result of our analysis is that 
(Vd, m /Vt,m) depends mainly on £d,o- I n the past, the con- 
clusions of various observational studies were that the disk 
and halo mass fractions in the optical parts of disk galax- 
ies depend on the SB and/or luminosity (or scale radius). 
Applying a direct linear multidimensional regression to 
(Vd, m /^4,m) with Ed,o and h d as independents, we find 
that: 



. Hi ) oc £ d ,o 



0.12±0.02t -0.07±0.03 



(7) 



with R=0.75. The large uncertainty in the slope of /id 
shows that the significance of this variable is small, in 
agreement with the result of the F-test. 

We arrive to similar conclusions by making use of 
the fundamental plane (FP) of disk galaxies. An exten- 
sive analysis of the luminous and baryonic FPs for our 
sample will be presented elsewhere. For the baryonic case 
and applying the direct linear regression fitting, we ob- 
tain that V t<m oc Ed ^^±0.02^0^0.03 ^ R=Q Q5 

Combining this relation with the definition of (Vd, m /Vt,m) 
(oc [Edfihd] 1 / 2 /Vt, m ), eq. is recuperated. 

We explore now any dependence of the (Vd, m /Vt,m) 
ratio on more direct observational quantities. We repeat 
the step-wise multiple regressions using now as indepen- 
dents LogEif.o and Loghx instead of LogSd,o and Logh d , 
besides (B — K) 1 T, and Log/ g . The statistical analysis 



shows again that only two independents are significant: 
Tik,o and (B — K) or Y,Kfi and T. The most significant 
correlation is: 



(Vd,m/Vt,m) OC T,k,0 



0.16±0.02 



(B-K) 



-0.08 ±0.02 



(8) 



with R=0.76. The regressions for Yik and T, and T,k 
and h K give: (V d , m /V t . m ) oc £ ^7±o.02 T o.oi4±o.oo4 
(R=0.75), and (%, m /V t , m ) oc E^o.nio.oi^-o.ofcto.oa 
(R=0.72), respectively. In the latter case, note the large 
uncertainty in the slope of h d . The trends seen above (eq. 
0) with the composite quantities Sd,o and h d remain ap- 
proximately similar when using the observables S^,o and 
h,K- We have also repeated the statistical analysis, this 
time neglecting the gas contribution, that is, Vd jm was cal- 
culated only from S s and h s (— h^). In this case, the cor- 
relation of (Vd, m /Vt,m) with the surface density becomes 
even slightly steeper than in the case including the disk 
gas parameters. Therefore, our conclusion about the de- 
pendence of (Vd,m/Vt,m) °n the disk surface density and 
its small dependence on the scale radius is robust w.r.t. 
the the assumptions made to calculate £d,o and h d from 
the observational data (see §2.2): the stellar M/L ratio 
7k, and the disk gas parameters, inferred from the total 
gas mass. 

It is important to remark that if the galaxy sample is 
restricted to a relatively small range of luminosities and/or 
SBs, then apparent anti-correlations between (Vd im /14 jm ) 
and h d , and luminosity or mass appear (in Figs. 7 and 8 
we show related effects for the mass-to-luminosity ratio). 
Therefore, for samples limited in luminosity and/or SB, 
the (Vd,m/14,m) ratio would seem to be larger for smaller 
galaxies. In fact, even for a complete sample, a slight trend 
in this direction appears, but it is probably due to the de- 
pendence of (Vd,m/14,m) on £d,o and the fact that smaller 
galaxies are typically also of lower SB. 

3.4. Observations and models: maximum disk? 

Observations and models show that the (Vd, m /Vt,m) ratio 
changes among galaxies mainly depending on their surface 
densities (brightnesses): it decreases monotonically as the 
galaxy is of lower surface density (Fig. 3a) . The horizontal 
lines in Fig. 3a show the range of maximum disk accord- 
ing to Sacket (1997), (V d /V t ) 2 2hd = 0.85 ± 0.10, which 
for the (Vd,m/Vt,m) ratio should be considered only as an 
approximation for low SBs). Notice that some fraction of 
the observed galaxies of ultra-high and high SB are de- 
scribed by maximum or marginally maximum disks, while 
LSB galaxies tend to be sub-maximum disks. Indeed sev- 
eral detailed observational studies of rotation curve shapes 
and SB profiles for HSB galaxies show that they tend to 
be maximum disks (e.g., Corsini et al. 1998; Salucci & 
Persic 1999; Palunas & Williams 2000). Theoretical ar- 
guments, including the swing amplifier constraints or bar 
slowing down due to dynamical friction also suggest a do- 
minion of luminous matter in the center of HSB galax- 
ies (see Bosma 2002 for a review and references therein). 
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On the other hand, for LSB galaxies, a maximum disk 
solution demands too high stellar M/L ratios, which for 
reasonable stellar population models imply non-realistic 
color indexes for the disks (de Blok et al. 1997,2001). The 
results in Fig. 3 complement all these previous works, uni- 
fying HSB and LSB galaxies in a continuity of decreasing 
luminous-to-total mass ratio within the optical part of the 
galaxy as the disk SB decreases. Nevertheless, this ratio 
varies significantly even for a fixed SB. Below we explore 
whether other galaxy properties play any significant role 
in defining the luminous-to-dark matter ratio. 

From the theoretical side, for disks formed within 
ACDM halos assuming /d= 0.05, only those with the 
highest surface densities could be described by the 
maximum disk solution. Galaxy models with moder- 
ate Ed.o (~ 400M Q pc~ 2 ) have sub-maximal disks, 
with (Vd,m/14, m ) « 0.66 ± 0.06 (for the NFW halos, 
(Vd, m /Vt,m) « v rf/ v t) 2 .2/i d a l most ever > as shown in Fig. 
2d). If /d is decreased to the more realistic values of 0.03- 
0.02 (Avila-Rccsc, Firmani & Zavala 2002), then models 
result with even lower (Vd, m /Vt, m ) ratios. The arrows in 
panel (b) of Fig. 3 indicate how much the (V^m/Vt.m) ra- 
tio of the models would decrease if instead of fd= 0.05, 
/d= 0.03 is used. In this panel the means per £ d ,o bins and 
their corresponding standard deviations are also plotted. 
The shaded region is for observations, while the region 
within solid lines is for the models. The thick dotted and 
solid lines connect the bin means of both groups, respec- 
tively. 

Observations and models show a large and roughly 
similar scatter in the (Vd, m /Vt im )-£d,o pl°t (Fig- 3). For 
the former, part of this scatter is due to the combina- 
tion of observational uncertainties in the measurements of 
V4, m , Ek,o, hx, and Mm propagated in the calculations 
of (Vd, m /Vt,m) and £d,o ( see Graham 2002 for uncertainty 
estimates for some of the observational data) . The resid- 
uals of the observational (Vd, m /Vt, m )-£d,o relation do not 
correlate in a significant way with any other galaxy pa- 
rameter, like / g , Mbar, Vt,m, morphological type, integral 
(B — K) color. Nevertheless, we find some trend for redder, 
earlier type, and more massive galaxies to have on average 
smaller (Vd, m /^t,m) ratios. The trends with the type and 
color would increase slightly if one takes into account that, 
as the galaxy is of earlier type and redder, the thicker its 
disk is. The thicker the disk, the lower the (Vd. m /Vt,m) 
ratio is. Notice that the observations have been uniformly 
corrected for disk thickness for all the galaxies (5% in 
(Vd,m/^t,m)). The trends just mentioned can be explained 
by our models. 

For the models, the scatter in the (Vd, m / 1 14 im )-Sd.o re- 
lation for /d constant is produced mainly by variations in 
the halo structure (concentration) and the mass (see Fig. 
2b). Since more concentrated halos are the product of ear- 
lier mass assembling (Avila- Reese et al. 1998), the disks 
formed within these halos with a gas accretion rate propor- 
tional to the MAH, will be redder. Thus, for a given Ed.o, 
redder disks are expected to have lower (Vd, m /Vt, m ) ratios. 



The weak correlation with the mass was already explained 
in §§3.1. These effects can be interpreted as weak devia- 
tions from homology and show a connection between halo 
structure and luminous evolution. Unfortunately, as men- 
tioned above, the observational uncertainty is too large to 
allow a direct comparison of models and observations re- 
garding the correlations of the residuals. Nevertheless, at 
least the trends seem to be in the same direction, showing 
the good predictive capability of the galaxy evolutionary 
models of FA2000. 

The comparison of observations and models in pan- 
els (a) and (b) of Fig. 3 shows some agreement among 
them, but on average models fall below observations. For 
fd = 0.03, the difference increases (see the arrows in Fig. 
3b). In panel (a) we also reproduce the line correspond- 
ing to the pseudo-isothermal model with /d= 0.05 (M v = 
3.5 1O 11 M , r c =3.9 kpc) presented in Fig. 3a, and add 
results for two other masses: M v = 3.5 10 10 M Q (r c =1.9 
kpc, upper line) and M v — 3.5 10 12 M Q (r c =7.7 kpc, lower 
line); the core radii are in accordance to inferences by 
Firmani et al. (2001), see §3.1. As was noted previously 
(Figs. 2a and 2b), the major difference between cuspy, 
gravitationally contracted, NFW and pseudo-isothermal 
halos in the (Vd !m /T4, m )-£d,o plane is at high surface den- 
sities. The observations lie actually in between both halo 
models, specially at high surface densities. 

The galaxies formed within ACDM halos seem to be 
somewhat less dominated by disk than observations sug- 
gest, specially when /d=0.03. The average difference in 
this case would be roughly of 10%. On the other hand, a 
pseudo-isothermal halo yields a (Vd, m /Vt, m ) ratio larger 
than observational inferences for HSB galaxies (and a 
steeper increase of the (Vd,m/Vt,m) ratio with £d,o)- These 
differences might decrease by using smaller core radii than 
assumed here or by introducing the disk gravitational 
drag over the pseudo-isothermal halo. The fact that the 
(Vd,m/14,m) ratio continues increasing with Ed,o for ultra- 
HSB observed galaxies could be pointing to that either 
the disk gravitational drag on the halo is not as efficient 
as in the models (the gravitational drag is the responsible 
that (Vd, m /Vt,m) tends to be constant as Sd,o increases to 
high values) or to that some physical process is expand- 
ing continuously the halo core (self-interacting CDM for 
instance). 

We conclude that the predictions of the ACDM sce- 
nario regarding the amounts of dark and luminous mat- 
ter in disk galaxies are in marginal agreement with ob- 
servations, specially if fd < 0.05. However, the way in 
which the amounts of dark and luminous matter change 
as different galaxy properties vary (Sd,o in particular) is 
in good agreement with observations. It seems that mi- 
nor modifications to the inner structure of the ACDM ha- 
los is enough to conciliate models with observations. In 
FA00 we have shown that the introduction of soft cores in 
agreement with inferences from the rotation curves of LSB 
galaxies, produce a ~ 10% increase of the (Vd,m/Vt,m) ra- 
tio (on average) for HSB galaxies. These models would be 
in better agreement with observations in Fig. 3. 
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4. Dynamical-to-baryonic mass and mass-to-light 
ratios 

In the previous Sect, we have investigated the (Vd, m /Vt,m) 
ratio from observed and model galaxies and determined 
how it depends on different galaxy properties. As shown 
in Fig. 2, panels c) and d), (Kl,m/^,m) seems to be good 
tracer of the amounts of luminous and dark matter in 
disk galaxies. In this section we will study directly the 
dynamical-to-baryonic and dynamical-to-stellar mass ra- 
tios, Mdyn/Mbar and Af dyn /Af s , as well as the mass-to- 
light ratios. The latter is the ratio most directly obtained 
from observations and it has been determined previously 
in the literature, mainly in the B— band (e.g., Salucci, 
Ashman & Persic 1991; Karachentsev 1991; Roberts & 
Haynes 1994; Graham 2002). We remark that Lk and 
Lb refer to the total galaxy luminosities, i.e. disk+bulge. 
Therefore, M s and Afb ar are the total disk+bulge stel- 
lar and baryonic masses, respectively. In general, while 
the (1^d, m /Vt,m) ratio refers to the disk contribution (Vd, m 
has been calculated from the disk SB profile (eq. [S]), the 
Af dyn /Afbar ratio refers to the total galaxy contribution 
(disk + bulge). Thus, from the observational point of view, 
the information contained in Af dyn /A/b a r i s slightly differ- 
ent to that contained in (Vd, m /Vt,m)- 

In §3, eq. Q}, we have shown the theoretical rela- 
tion between our observational quantity (Vd, m /Vt,m) and 
Mdyn/Mh&r defined at r pca k, the radius where the rota- 
tion curve peaks. For disk dominated galaxies, the radii 
where the disk and total velocities peak are approximately 
the same, at 2.2ft, d (% = 2.2), and the disk mass fraction 
at this radius is Jl = 0.64. For halo dominated galax- 
ies, the departure of x and Jl from 2.2 and 0.64 will de- 
pend mainly on the halo mass distribution; in any case, 
both will increase, in such a way that the ratio xj /z, will 
only slightly change (increase). For example, for NFW ha- 
los, x//l remains almost constant (« 2.2/0.64) for al- 
most any case (see Fig. 2d). Thus, we may estimate ap- 
proximately the Mdyn/-Mbar ratio at 2.2/i d by multyplying 
M dya /M bax = (V t ^ 2 h d /G)/M har by 2.2/0.64=3.44. We 
may also estimate Af d y n /Afb a r at 5/i d = 5r pGa k/x (roughly 
the optical radius of disks). For this, we have to assume 
that the rotation curve remains flat from its maximum up 
to this radius, and M dyn should be multiplied by 5; for an 
exponential disk, the mass at 5/i d is close to the total disk 
mass (99%). In fact, the outer slope of the rotation curves 
changes, mainly depending on the disk surface density and 
luminosity. While some LSB and low-luminosity galaxies 
may have still increasing rotation curves at ~ 5/i d , HSB 
and luminous galaxies may present decreasing shapes at 
this radius (e.g., Casertano & van Gorkom 1991; Persic et 
al. 1996; Verheijen 1997). Galaxies with increasing rota- 
tion curves at the last measured point were excluded from 
our sample. The highest SB and luminosity galaxies may 
have on average a rotation velocity at 5/i d smaller than 
V t , m by a 10 - 20%. Therefore, the M dyn /Af ba r ratio at 
5/id presented below could be overestimated by this factor 
for the highest SB and luminosity galaxies. 



4.1. Results 

In panels (a), (b), and (c) of Fig. 4 we plot 
(Vt, m 2 h d /G)/M haI vs. E d , , {V t , m 2 h d /G)/M s vs. E s>0 , and 
(Vt, m 2 hd/G)/LK vs E^o, respectively. The axis to the 
right of each panel was shifted up by a factor of 5 to resem- 
ble approximately M dyn /M b ar, Af dyn /M s , and M dyn /LK 
at 5 hd, respectively. To obtain these ratios at 2.2 hd, the 
shift should be approximately by a factor of 3.44. Similar 
to (Vd, m /Vt im ), the Md yll /-Mbar ratio presents a depen- 
dence on Ed,0i with a poor correlation with ft, d . Applying 
a direct multidimensional linear regression to Af dyn /Afbar 
as a function of Ed,o and hd (Log-Log), we arrive to the 
same conclusion of the previous section: the dynamical- 
to-baryonic mass ratio in disk galaxies depends mainly on 
£d,o but weakly on scale (hd or Afbar/- The correlation 
obtained from the observations is 



M dyn /M h 

ill oc E d n 



-0.31±0.03i -0.05±0.05 



(9) 



with R=0.8. This is consistent with eq. |Q if Af dyn /Mb a r oc 

(v d , m /v t , m )- 2 . 

The direct linear regression of Log(Afd yn /Mb ar ) with 
LogEd.o gives M dyn /M bar = al0 - 63±010 E d , -°- 31±0 04 
(R=0.74), with a = 5 (solid line in Fig. 4) and a = 2.2 
for this ratio defined at 5 and 2.2 hd, respectively. The 
ratio of baryon to dark matter, Afb ar /Af dm , is related 
to A/d y n/A/bar by: Af bar /M dm = l/(Af dyn /M bar - 1). 
Thus, on average the (Afb a r/Afdm) 2 2/i d ra ti° f° r ultra- 
HSB (E d ,o ~ 2OOOM pc- 2 ), moderate-HSB (Z d ,o « 
4OOM pc- 2 ), and very-LSB (E d: o ~ 5OM pc~ 2 ) galaxies 
are 2.5, 0.8 and 0.25, respectively. For the MW and M31 
galaxies, (Afbar/Afdm)2.2 = 0.74 and 1.34, respectively. 

The anti-correlation of Af dyn /Mb a r with the surface 
density becomes steeper (and tighter) when one passes 
from the baryonic to the stellar disk or the K— band 
(panels b and c). A direct multidimensional lineal regres- 
sion shows that M dya /M s oc E s Q -o.5i±o.o4^-o.i3±o.06 
and M dyn /L K oc Sjfj() -o-B6±o.03^-o.X6±o.06 ( R=0 , 86 and 
0.89, respectively). This change is because galaxies with 
lower surface densities have typically higher gas fractions 
than the higher surface density ones. A contrary effect 
happens when one passes from stellar masses and densities 
to the B— band luminosities and SBs: lower surface den- 
sity galaxies are typically more luminous in the B— band 
(bluer) than the higher surface density galaxies with the 
same masses. Therefore, the Md yn /LB vs. E_b,o relation 
(Fig. 5a) becomes shallower and noisier than the Af dyn /M s 
vs. E s .o relation. Nevertheless, we find that Af dyn /is still 
(anti) correlates significantly with T<b,o, and does not cor- 
relate with hs (Fig. 6a). A multidimensional linear regres- 
sion shows that 



i\/f it ^ V -0.37±0.04, 0.03±0.06 
M d yn/^B OC L,Bfi tlB 



(10) 



with R=0.77. 

It should be noticed that if the sample is strongly re- 
stricted in magnitudes, a spurious correlation of Af dyn /is 
with hs appears. For instance, if the data in Fig. 6a 
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Fig. 4. (a). Dynamical-to-baryonic mass ratio vs. the 
baryonic surface density for observations (empty circles) 
and models (same symbols as in Fig. 3). The left axis is 
the quantity inferred directly from the observational data, 
while the right axis is this quantity multiplied by 5; in the 
assumption that the rotation curve remains constant up 
to 5/id, this corresponds to M dyn /M bar at 5h d . Solid and 
dotted lines are the linear regressions for models and ob- 
servations, respectively. For observations, the fitting (at 
5/ia) is showed in the upper right, (b). Same as in panel 
(a) but for the stellar mass and surface density, (c). Same 
as in panel (a) but for the if— band luminosity and SB. 
Only observations are showed in this panel. 



(Mdyn/LB vs. hs) are sorted into three groups according 
to their B— band magnitudes, then each group show an 
apparent correlation. We also have calculated the residu- 
als of the M<± yn / Lb-^b,o and /ib-£b,o relations. There is 
not any dependence among them. Contrary to this, and 
using the FP, Graham (2002) arrived to the conclusion 



that Mdy n /LB does not depend on ~Eb,o but with Kb- In 
§4.2 we discuss the reasons of this apparent difference. 

The correlations we have found for M dyn /LK and 
T,k,o, and for M^u/Lb and £b,o imply the existence of a 
dependence of L^/L^ on both Eb,o/^k,o and T,k,o- With 
some scatter, this dependence can be written as 



(B-K)?a 0.36(B - K) - 0.19/^,0, 



(11) 



which agrees with a direct comparison of these observable 
parameters. The disk center is typically redder than the 
average disk color and the difference increases for redder 
galaxies (e.g., de Jong 1996c; Vcrhcijcn 1997). 

In Fig. 7a we explore whether the M dyn /M har ratio 
(defined at 5h d ) correlates with the integral (B — K) 
color. Apparently there is no significant direct correla- 
tion, but a weak anti-correlation (see also Graham 2002). 
This is mostly because redder galaxies are commonly of 
higher SB. As was seen above, galaxies with higher SB are 
less dark matter dominated than galaxies with lower SB. 
However, for a given SB, one expects that redder galaxies 
(assembled earlier) are more dark matter dominated, at a 
given SB, than the bluer galaxies (Avila- Reese & Firmani 
2000). This is because dark halos assembled by an early 
collapse are more concentrated than halos assembled by 
extended accretion. In Fig. 7b, results from the same evo- 
lutionary models presented in §3 are shown in terms of 
the (B — V) color. The data were sorted into three groups 
according to their disk surface densities. For the lower sur- 
face density models, M dyn /Mb al correlates with the color 
as expected. For higher surface densities the correlation 
tends to disappear because of the effect of the surface den- 
sity and the "saturation" of the dark-to-baryonic mass ra- 
tio at high SBs. This is due to the disk gravitational drag 
over the halos (otherwise, the bluer HSB models could 
have lower M dyn /Mb al ratios). In panel (a) the observa- 
tional data has also been sorted into the same three groups 
according to their disk surface densities as in panel (b). 
The trend for the low surface density galaxies is similar 
than that for the models. At high surface densities obser- 
vations also agree roughly with models. However, several 
of the highest surface density and reddest galaxies appear 
with lower Md yn /Afb ar ratios than models. These galax- 
ies, having the most dominant disks and (probably) the 
more concentrated halos, reveal that dark matter is not 
as dominant as in the case of disks formed within NFW 
halos. 

In Fig. 4, panels (a) and (b), the data from the evolu- 
tionary models are also included (solid symbols; but notice 
that there are no model predictions in the K band). Solid 
and dashed lines are direct linear regressions to the ob- 
servational and model data, respectively. As in the case 
of the (Vd,m/Vt.m) ratio, models agree marginally with 
observations, yielding on average higher -Md yn /Mb a r and 
Md yn /M s ratios than observations, specially for the high- 
est surface density galaxies. For example, for the groups of 
ultra-HSB and moderate-HSB galaxies mentioned above, 
the model {M\ }av /M dm ) 2 ratios are on average 1.4 and 
0.6 (with a large scatter), which means approximately 1.8 
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Fig. 5. (a). Same as Fig. 4, but for B— band luminosity 
and SB, and only for observational data. (b). Same as 
panel (a) but using the definitions of Graham (2002; see 
text) and including some LSB galaxies with large bulges 
also considered by Graham (circles with crosses). Notice 
how the scatter in the correlation seen in panel (a) in- 
creases in panel (b). 



and f .25 times more dark matter dominated than obser- 
vations, respectively. For the LSB galaxies, the ratios are 
roughly similar. The differences between observations and 
models increase if /a is fixed to values smaller than 0.05. 
The halos of observed galaxies have on average less dark 
matter content in the center than the ACDM model sug- 
gests. It should be also noted that the differences between 
models and observations decrease when passing from stel- 
lar to baryonic (stars+gas) parameters. This is because 
models have slightly higher gas fractions than observa- 
tions (FAOO). In our models all the gas is cold, while the 
observations suggest that some fraction of the gas may 
be ionized in the outer disks. Therefore, this gas in not 
detected in the HI observations. 

4.2. Comparison with previous works 

Using the tilt of the optimal fundamental plane (FP) of 
disk galaxies, Graham (2002) concluded that Md yn /Ls 
depends on hs rather than on Hb,o- Since we share al- 
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Fig. 6. (a). Same as in Fig. 5a but for hs instead of 
Sb,o in the x-axis. The data are sorted into three groups 
according to their B-band magnitudes. The correspond- 
ing ranges are: 4.0 I0 8 — 2.5 10 9 Lb q (open symbols), 
2.5 I0 9 - 1.6 1O 1O L B0 (targeted symbols), and 1.6 10 10 - 
1.0 10 11 Lbq (solid symbols). Although for all the data is 
not significant, for the data sorted in luminosity ranges, 
a correlation between the mass-to-light ratio and Kb ap- 
pears, (b). Same as in Fig. 5b but for hg instead of Sb,o 
in the x-axis. A correlation appears for all the data in this 
panel (calculated according to Graham 2002). 



most the same observational sources, our results broadly 
agree with those of Graham (2002) if wc follow his FP 
analysis. In this indirect way to estimate the dependence 
of Mdyn/^s on Sb,o and hs, the multidimensional re- 
gression is applied not to Mdyn/Le, but to 14 ,m (the FP). 
This is one of the reasons behind the different conclu- 
sions achieved in Graham's paper. From his data, one 
sees that if the mass-to-light ratio dependence on Sb,o 
and Hb is determined directly ("kinematic" M/L ratio) 
instead of through the FP (the "photometric" M/L ra- 
tio), then the dependences on T,b,o and hs increases and 
decreases slightly, respectively (see his Fig. 4), and the re- 
sult approaches to what we have reported in the previous 
subsection. 

Other reasons to explain the difference in results are 
the way to calculate Md yn and Lb- Graham (2002) uses 
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the disk scale radius in the current band while we define a 
unique physical Md yn , using the baryonic disk scale radius 
/id- Rather than the measured Lb, Graham uses Sb^/ib 2 . 
Figures 5 and 6, panel (b), show our data using the same 
definitions as in Graham (2002). We also include in these 
plots LSB galaxies with large bulges from the sample of 
Beijersbergen (1999) and de Blok et al. (1998) (circles with 
crosses; for these galaxies there is available information 
only in the B— band), which were used by Graham. It is 
clearly seen that the dependence of M dyn / Lb on £_b,o be- 
comes weaker in panel (b) of Fig. 5 than in panel (a). On 
the other hand, as in Graham (2002), some correlation 
with hs appears now (compare panels a and b of Fig. 
6). Thus, our results agree broadly with those of Graham 
(2002) when we apply to the data his technique and def- 
initions, and when we include the peculiar sub-sample of 
LSB galaxies with large bulges. 

Nevertheless, we think that our analysis and the con- 
clusion in §4.1 that M dyn /LB anti-correlates with T.b,o 
and does not depend on Kb or Lb is physically correct. If 
we use the result of Graham (2002), it is easy to show that 
(B -R) « -0.25(5 - R) - (0.09^b + 0.12Mb), which 
implies an anti-correlation of the central color with the 
integral color of galaxies, contrary to that inferred from 
direct observational comparisons. 

For a sample of 114 spirals, Karachentsev (1991) has 
found that the ratio of dynamical mass-to-L b does not de- 
pend on \jb but does anti-correlate with the B— band SB, 
in agreement with our result (see also Giovanardi & Hunt 
1988). The anti-correlation in the data of Karachentsev 
(M dyn /L B oc S^j S q, s ~ 0.5) is steeper than we have found. 
We expect that the slope s would decrease if the velocity- 
dependent internal extinction correction applied here is 
applied to those data. 

Salucci et al. (1991) inferred the ratio M dyn /Mb al , by 
using B— band data for a sample of HSB galaxies and as- 
suming constant central SB. They concluded that this ra- 
tio decreases roughly as L^ 0A . However, as we have men- 
tioned in §3.3, a spurious correlation among (Vd., ro /Vt,m) 
and luminosity (mass) appears if the sample is restricted 
in SB. The same applies for the M dyn / ratio. Salucci 
et al. (1991) assumed constant SB and their sample was 
indeed limited to HSB galaxies only. Furthermore, we ex- 
pect that the correlation with Lb in Salucci et al. will 
diminish if a luminosity-dependent internal extinction co- 
efficient is applied to the data. 

Finally, from the analysis of a sample of optical rota- 
tion curves from HSB late-type spirals, Persic & Salucci 
(1988, 1990) have found that M dyn /M ba r at 3.2h d de- 
creases with Lb as L~^ A . The method is conceptually 
different from that used here and it is based on a different 
observational information. Our results (observational and 

theoretical) show that (-Md yn /Mbar) 2 2 ft, d ^ oes n0 * corre ~ 
late with Mb ar or Lb- The main assumption we have done 
to infer the M dyn /M\ )ar ratio at 2.2h d from observations 
is that x/ Jl{x) in eq. |Q| remains almost constant with 
the SB or luminosity. This is indeed true for disks formed 
inside CDM halos (or for any halo with a circular velocity 
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Fig. 7. (a). Approximate dynamical-to-baryonic mass ra- 
tio at 5/id vs. integral (B — K) color for observational 
data. No correlation, if any, an anti-correlation, is seen. 
However, for limited ranges of surface density, a corre- 
lation appears as this density is lower. The data in the 
plot were sorted into three groups according to their disk 
central surface densities: open, targeted and solid sym- 
bols are for low, intermedium and high surface densities, 
respectively (we use surface densities instead of Lk lumi- 
nosity in order to have the same ranges for the model data 
presented in panel [b]). The solid line is an eye-fit to the 
low surface density group, (b). Same as in panel (a) but 
for the model data and for (B — V) instead of (B — K). 
See text for explanations of this Figure. 

that does not increases significantly after the maximum of 
the disk component), where (Vd, m /Vt, m ) « (Vd/Vt) 2 2 f ld as 
seen in Fig. 2d. Since our analysis in Sect. 3 of the trends 
of (Vd,m/Vt,m) ° n several galaxy parameters has shown a 
reasonable agreement between observations and models, 
we may speculate that the galaxy halos do not depart sig- 
nificantly from the NFW mass distribution, in particular 
at intermedium radii. Nevertheless, more tests are needed 
in order to see whether x/Jl{x) indeed does not depend 
on SB/luminosity. 

A more direct comparison of the results from our ob- 
servational sample and the inferences of Persic and col- 
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laborators, can be done through the (Vd, m /Vt, m ) ratio, 
a quantity that was calculated from the available data 
without introducing any assumption about x/Jl(x). The 
analysis presented in §3.3 shows that (Vd, m /Vt, m ) for our 
sample does not depend significantly on mass or luminos- 
ity, but on SB. In the opposite, the mass decomposition of 
the universal rotation curve implies a strong correlation 
of (Vd,m/Vt, m ) with luminosity (see Fig. 6 in Persic et al. 
1996), while the correlation with SB was not explored in 
that work. 

The investigation in detail of the differences on the re- 
sults obtained with the two different methods is out of the 
scope of the present paper. Our results suggest that the SB 
should be taken into account in the kinematical method 
as an important parameter for the shape of the rotation 
curves. On the other hand, we notice that, according to 
our sample data, the /is //id (or hs/hic) ratio increases on 
average with Lb- The results of Persic & Salucci (1988, 
1990) and Persic et al. (1996) are based on the determi- 
nation of the slope of the rotation curve at 3.2 scale radii. 
They used the optical (mostly B— band) disk scale radii 
instead of ha or h^; the latter are more directly related 
to the gravitational field of the disk than hs- Since the 
slope of the rotation curve it is sensible to the radius at 
which is measured, a systematic shift of this radius toward 
larger values as the luminosity increases, will introduce an 
(extra) anti-correlation of this slope with the luminosity, 
if the the rotation curves are actually more decreasing as 
the galaxy is more luminous or of higher SB. We specu- 
late that, by using the more physical Hk or radii rather 
than hs or hj, the anti-correlation of the rotation curve 
slope at 3.2 scale radii with Lb found by the mentioned 
authors will become weaker. 

In general, most of the previous works on the dark and 
luminous contents in disk galaxies did not take into ac- 
count that these contents vary significantly among galax- 
ies, depending on the disk SB. We have shown that the 
ratio of total-to-luminous matter content within the opti- 
cal regions of the galaxies decreases as the SB increases, 
without any significant (anti)correlation with galaxy scale. 

5. Summarizing discussion 

We have presented a compilation of normal disk galaxies 
with available photometric parameters in both the B and 
K bands, rotation curves and/or HI- velocity line- widths, 
and HI integral fluxes. A total of 78 galaxies have been se- 
lected as useful for our analysis. We have uniformized the 
data (luminosities, disk central SBs and scale radii, maxi- 
mum rotation velocity, disk gaseous mass) and used them 
to further calculate several stellar and baryonic galaxy pa- 
rameters. With this sample we explored the ratios of dark 
to luminous mass (and luminosity) within disk galaxies 
and how these ratios vary with the galaxy properties. The 
inferences obtained here are alternative to other more di- 
rect methods, as the rotation curve decomposition and the 
statistical analysis of rotation curve shapes. A detailed 
comparison of our results with those of more direct meth- 



ods is not within the scope of the present paper, but de- 
serves a future analysis. In this paper, our observational 
results were compared with model predictions and ana- 
lyzed further in the light of the behavior of the models. 
We next discuss the main conclusions of this work. 

• The (Vd, m /Vt jm ) ratio inferred from observations de- 
pends mainly on the disk central surface density, Xd.o (or 
T,k,o)i increasing approximately as £d,o ' 15 - Only a frac- 
tion of HSB and ultra-HSB galaxies in our sample fall in 
the range of maximum disk ((Vd, m /Vt, m ) ~ 0.85 ± 0.10); 
as the disk central surface density decreases, the galax- 
ies become more dark matter dominated, with the low- 
est SB galaxies having (Vd, m /Vt, m ) ratios around 0.5. The 
MW and M31 galaxies have {Vd,m/V t , m )= 0.67 and 0.84, 
respectively. In previous works, by using the analysis of 
rotation curve shapes, a "democratic" solution to the dis- 
cussion about maximum or sub-maximum disk was found, 
in the sense that there are galaxies with rotation curves 
highly dominated by the disk component and other ones 
dominated by the halo component (e.g., Persic et al. 1996). 
The main driving parameter in these works was the lumi- 
nosity or mass. Our results show instead that the driving 
parameter is the disk surface brightness or density. 

• Although the (Vd, m /Vt im ) ratio depends mainly on 
Sd,o (^k,o), a large scatter is observed in the (Vd, m /Vt, m )- 
^d,o plot- Interestingly, we find some trend of decreas- 
ing (Vd,m/Vt, m ) as the galaxy is redder, of earlier type 
and more luminous (massive) . These results hint us about 
a connexion between the luminous and morphological 
galaxy properties to the dark halo properties. The latter 
ones are related to the initial (cosmological) conditions. 

• To gain a deeper understanding of the observations, 
we have confronted them with models of disk galaxy for- 
mation and evolution in the ACDM cosmogony (FA00). 
For a disk-to-total halo mass ratio /a= 0.05, the mod- 
els agree marginally with the observational data in the 
(Vd,m/Vt, m ) vs. Sd,o plane (Fig. 3). Models have on aver- 
age slightly lower (Vd, m /Vt, m ) ratios than observations, the 
difference becoming larger for the highest values of Sd,o- 
For /d=0.03, the differences increase slightly. These differ- 
ences suggest that real galaxies may have less dark mat- 
ter in the center than disks formed within ACDM halos. 
The observations in the (Vd, m /Vt, m ) — Sd,o plot lie indeed 
in between models with cuspy NFW halos and pseudo- 
isothermal halos with shallow cores, specially at the high- 
est surface densities, where the differences between both 
halo models become large. Nevertheless, the differences 
between observations and the ACDM models are modest, 
and the dependences of (Vd, m /Vt, m ) on different galaxy 
properties (in particular Sd,o) for observations and these 
models agree rather well (e.g., see thick dotted and solid 
lines in Fig. 3b). Therefore, only a minor modification to 
the inner structure of the ACDM halos seems to be nec- 
essary to remedy the differences with the observations. 

A possible mechanism to form soft cores in the CDM 
halos, without altering their MAHs and outer struc- 
tures, is the gravothermal expansion produced in the 
core if the CDM particles are allowed to self-interact 
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among them. Cosmological N-body simulations for self- 
interacting CDM particles with a modest cross section in- 
versely proportional to the particle relative velocity, have 
shown to be an attractive alternative to produce soft cores 
in agreement with observational inferences from dwarfs to 
galaxy cluster scales, keeping at the same time the success- 
ful predictions of the "standard" ACDM paradigm (Colin 
et al. 2002). The constant expansion of the core could ex- 
plain why the (Vd, m /Vt,m) ratios of ultra-HSB galaxies is 
so high in spite of the strong gravitational drag of the disk 
on the halo (see Fig. 2). Models of disk galaxy evolution 
within the self-interacting ACDM halos will show whether 
self-interaction leads to a better agreement with observa- 
tions than in the case of collisionless CDM particles. 

• The models show that, in spite that more massive (lu- 
minous) systems have typically less concentrated dark ha- 
los than the smaller ones, for a given Sd,o, the former are 
more dark matter dominated (smaller (V^m/Vt.m) ratios) 
than the latter. This trend, also seen in the observations, 
is easily explained in the light of the models (§§3.1, 3.2.2). 
On the other hand, for a given mass, the (V^m/V^m) ra- 
tio of disks in more concentrated halos is lower than in 
the less concentrated ones. More concentrated halos are 
typically the result of earlier mass assembling, therefore 
the disks formed within them are redder because the gas 
infall rate at latter epochs is low. Thus, the trend of red- 
der galaxies presenting smaller (Vd, m /Vt,m) ratios, at least 
for LSB galaxies, is understandable in the models and it 
might explain the same trend found for the observations. 
In general, the scatter of the observational data in the 
(Vd, m /^t,m)-Sd,o plane is only slightly larger than for the 
models (Fig. 3b). Most of this scatter is explained in the 
light of the models, and the rest is perhaps due to the 
observational uncertainties. 

• We also have presented the dynamical-to-baryon 
mass ratio of models and observations (Fig. 4). To esti- 
mate this ratio at 2.2 and 5 x /id from the kind of obser- 
vations we have used here, one has to assume that Vt, m 
occurs at 2.2 hd for all the galaxies, or that 14, m remains 
constant up to bha, respectively. The Md yn /Mb a r ratio of 
observations and models depends mainly on Ed.o an d does 
not depend on scale (hd or A/bar) ■ These trends sharpen 
when using (M s , £ s ,o) and (Lk, ^k,o) instead of (Mb ar , 
£d,o)- But when passing to (Lb, Sgj), the relation be- 
tween Md y n/L>B and £_b ; o becomes shallower and noisier 
than in the other cases (Md yn /M s vs. £ s ,0j for example). 
Nevertheless, we find that Md yn /£s still (anti)correlates 
significantly with ~Sb,o, and does not correlate with Kb 
(or Lb)- The Md yn -to-Mb a r, -M s , -Lk, and -Lb ratios at 
5/id for the lowest SB galaxies are approximately 8, 14, 
10.5 and 18, respectively, while for the highest SB galax- 
ies are 2, 2, 1.5, and 5, respectively (Figs. 5 and 6). The 
Mdyn/Mbai ratio for the highest SB model galaxies (max- 
imum difference with observations) is around 2.4. This 
means that for the highest SB galaxies, models predict 
halos ~ 1.7 more massive at 5hd than inferred from ob- 
servations (at 2.2 hd this difference is of ~ 1.8). 



• The Mdyn/Mbar ratio does not correlate with Lb (or 
Mb ar ) as was found previously in works based on the sta- 
tistical analysis of rotation curve shapes, where the deter- 
mination of the rotation curve slope at some typical radius 
(3.2hs in these works) is crucial. We have noted that the 
correlation of hs/hd (or hs/ha) on Lb could contribute 
significantly to an anti-correlation of the slope with Lb, 
in the understanding that hd (or hx ) is the physical ra- 
dius to be used. We also note that the strong dependence 
of the (Vd,m/T4,m) ratio on luminosity inferred from the 
universal rotation curve is not reproduced by our data 
analysis. The point is that two different methods, using 
different samples, arrive to different conclusions: we find 
that the the fractions of dark and luminous matter within 
the disks depend mainly on the SB or surface density, 
while in works based on the analysis of the rotation curve 
shapes, the main parameter is the luminosity or mass. 
In fact luminosity and SB are correlated, but the deter- 
mination of which of them is the dominant is an impor- 
tant issue to understand galaxy formation and evolution. 
If the luminosity (or mass) is the main parameter, then 
disk formation in ACDM models could not be realistic or 
intermediate astrophysical processes like a strongly mass- 
dependent feedback could be relevant. It is important that 
future works solve the controversy stated here. 

• The correlations of the mass-to-light ratios with SB 
in the B and K bands found here imply that the inte- 
gral color (B — K) correlates with the central disk color 
(B—K)o fea.|llp. in accord to a direct comparison of both 
parameters. We have also found that as galaxies are of less 
SB, a more significant correlation among the Mdyn/Mbar 
ratio and (B — K) appears (Fig. 7). According to the mod- 
els, galaxies are redder in part because they formed in 
more concentrated dark halos (halos that collapsed earlier 
on average). Therefore, for a given SB, Md yn /Mb a r tends 
to be larger (more dominion of dark matter) for redder 
galaxies. However, for disks with high SBs, the Md yn /Mb ar 
ratio tends to be more constant than for low SBs (due to 
the disk gravitational drag over the halo), in such a way 
that the correlation with the color minimizes. The similar 
features found between models and observations (compare 
panels a and b of Fig. 7) suggests that this is a plausible 
explanation for the observational trend. 

• The good agreement found between models and ob- 
servations in the comparisons presented along the present 
paper favors the ACDM scenario, under the assumptions 
that disk galaxies are formed by gas accretion rather than 
mergers, and that during the collapse of the gas, detailed 
angular momentum conservation is obeyed. In the light 
of the results presented here, the main difficulty of this 
scenario is an apparent excess of dark matter within the 
disk optical radius. A minor modification to the physics 
of the formation of the ACDM halos seems to be enough 
to overcome this difficulty. 
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Appendix A 

In this Appendix we describe the construction of 
the disk+halo composite models. The overall procedure 
closely resembles that one presented by Mo et al. (1998). 

The disk is assumed to be in centrifugal equilibrium, 
to be thin, and to have an exponential surface density 
distribution, E^oM = ^d.o^xp{—r / h d ) , where the central 
surface density £d,o and scale radius h d are related to 
the total disk mass by Mb ar = 2irT, d _oh d 2 . The circular 
velocity is calculated according to eq. ©. The angular 
momentum of the disk, neglecting its gravitational effect, 
is: 



J d = 2tt / V h (r)^ Q (r)r 2 dr, 



(12) 



where V/j(r) is the halo circular velocity. It is assumed 
that the disk forms from the baryons trapped within the 
virialized halo (with a mass fraction /d) and that they 
follow originally the same mass and angular momentum 
distribution of the halo. If detailed angular momentum 
conservation is obeyed, then J d = J, where J is the halo 
angular momentum, and it is related to the dimensionless 
spin parameter A by: 



(13) 



where E is the total energy of the halo. Now, supplanting 
these relations into eq. (|12fl . the disk scale radius h d is 
calculated; it will be proportional to A. Therefore, given 
M v , /d, A and the halo density profile, the two parameters 
of the exponential disk in centrifugal equilibrium remain 
defined. Following, two halo models are analyzed: 

Pseudo-isothermal (P-I) sphere. The density and cir- 
cular velocity profiles of a P-I sphere with core radius r c 
and constant central density po are given by: 



p{r) 



Po 



1 + (r/r c ) 2 



= [4ttGVc(1 - — arctan(— ))] (14) 
7* r c 



The asymptotic velocity is V aS ym = {^Gp^r 2 ) 05 . 
According to the cosmological spherical collapse model, 
at the epoch of virialization of the last shell of an over- 
dense homogeneous sphere, its radius (virial radius, r v ) 
is such that < p(< r v ) >= A c p u , where < p(< r) > 
is the average density of the sphere at r, A c is a con- 
stant, and p u is the matter density of the universe at that 
epoch. For the present epoch and for the ACDM cosmol- 
ogy used here, A c = 334 (Eke, Navarro & Frenk 1998). 
The mass contained within r v is the virial (halo) mass, 
M v . Taking into account that r„/r c >> 1, it is easy to find 
that V asym « {GM v /r v )^ 2 « 0.014(Af v /M Q ) 1 /3k mS - 1 , 
and that E k, — (M v V 2 sym ) /2. At radii larger than r c , 



Vh rapidly tends to V asym =const, and the integral in eq. 
(I12f> can be approximated by Jd w 2Mb ar /idV asym . Now, 
from Jd = J and using eq. (|13fl . ^d can be calculated, 
h d x \M v 1/3 f d -\ and S d! o oc M v 1/3 f d 3 \~ 2 . 

Once M v , /d, and A are defined, we calculate the 
disk parameters h d and Sd,o and use them in eq. JSJ 
to get the disk velocity component. The halo velocity 
component is calculated according to eq. @2J. The core 
radius r c is taken from empirical inferences found by 
Firmani et al. (2001) from dwarf and LSB galaxies: r' c ~ 
5.5(V asym /100kms -1 ) - 95 , where we have to take into ac- 
count that the core radius of the non-singular isothermal 
model used by them, r' c , is 1.45 times the core radius of 
the pseudo-isothermal model used here, r c = r' c /1.45. 

NFW cosmological halo model. The typical density and 
circular velocity profiles of the CDM halos are given by: 



p{r) 



(r/r s )(l + r/r s 



GM h (r) ^Gp s r 3 [fn(l + £) - 



(15) 



where the parameters p s and r s are related. Defining 
again as above the virial radius r v and introducing the 
concentration parameter c = r v /r Sl one obtains that 
p s = {A c p u /3)c 3 g(c), where g{c) = l/(/n(l+c)-c/(l+c)). 
Now, the halo circular velocity can be expressed in terms 
of the velocity at the virial radius (V v =(GM v /t v ) 1 ' 2 ): 
(Vh/Vv) 2 = (l/w)(ln(l + cw) - cw/(l + cw))g(c), where 
w = r/r v . As shown above, from the spherical collapse 
model one obtains that V v oc M^ 1 ^ . Therefore, given M v 
and c, the NFW velocity profile remains defined. In or- 
der to calculate the disk parameters, the total halo en- 
ergy is necessary: E = —GM v 2 f(c)/(2r v ), where /(c) is 
a function of the concentration parameter (see Mo et al. 
1998) and r^ ocV„. Cosmological numerical simulations 
show that c slightly decreases with M v and has a signifi- 
cant scatter (AlogcNFW ~ 0.18, Bullock et al. 2001). 

Gravitational contraction of the NFW halo. The for- 
mation of a dense disk inside the cosmological dark halo 
induces a gravitational contraction of the halo in the inner 
regions. Under the assumption of radial adiabatic invari- 
ance, one may calculate the mean radius r where ends 
a particle with an initial mean radius r.^: GMf(r)r = 
GM(n)n, where M f {r) = M d {r) + M{n)(l - f d ) is 
the final total mass distribution and M(ri) is the ini- 
tial halo (NFW) mass distribution. To calculate the disk 
mass distribution Md(r), h d should be found as described 
above. But to calculate h d the total circular velocity 
V 2 (r) =V^(r)+V 2 ! (r) should be given, where Vj t (r) = 
G[Mf(r) — Md(r)]/r. Given M v , f d , A and c, an iterative 
procedure yields h d and Vt(r) after the halo contraction. 
This procedure is explained in Mo et. al (1998) and we 
refer the reader to that paper for more details. 

Appendix B. 

Here we present the main ingredients of the self- 
consistent evolutionary models to be studied in this pa- 
per. The disk is built up within a growing ACDM halo. A 
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special extended Press-Schechter approach is used to gen- 
erate the statistical mass aggregation histories (MAHs) of 
the halos from the primordial density fluctuation field, 
and a generalized secondary infall model is applied to 
calculate the virialization of the accreting mass shells. 
The evolution and structure of the ACDM halos calcu- 
lated this way agree well with results from cosmological 
N-body simulations (Avila- Reese et al. 1999); halos assem- 
bled through early active MAHs end more concentrated on 
average than halos with extended MAHs. The mass shells 
are assumed to have aligned rotation axis with specific 
angular momentum given by j s h(t v ) = dJ(t v )/dM v (t v ), 
where J = \GMl /2 / \E\ l/2 , J, M„ and E are the total 
angular momentum, mass and energy of the halo at the 
shell virialization time t v . The spin parameter, A, is as- 
sumed to be constant in time. As the result of the assem- 
bling of these mass shells, a present day halo ends with an 
angular momentum distribution close to the (universal) 
distribution measured by Bullock et al. (2001) in N-body 
simulations. A fraction /d of the mass of each shell is as- 
sumed to cool down and form a disk layer in a dynamical 
time. The radial mass distribution of the layer is calcu- 
lated by equating its specific angular momentum to that 
of its final circular orbit in centrifugal equilibrium (de- 
tailed angular momentum conservation is assumed). The 
suppcrposition of these layers form the disk. The gravi- 
tational interaction of disk and halo is calculated using 
the adiabatic invariance formalism. The local SF is trig- 
gered by the Toomre gas gravitational instability criterion 
and self-regulated by a vertical disk balance between the 
energy input due to SNe and the turbulent energy dis- 
sipation in the ISM. The SF efficiency depends on the 
gas surface density determined mainly by A, and on the 
gas accretion rate determined by the cosmological MAH. 
Finally, we calculate the formation of a secular bulge by 
using the Toomre instability criterion for the stellar disk. 
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Table 1. Raw and corrected photometric and dynamical parameters of the disk galaxy sample 
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a Luminous distance in Mpc. 

Raw uncorrected apparent magnitude. 
c Raw uncorrected disk central SB (magarcsec -2 ). 
d Corrected integral (B — K) color. 
c Disk scale radius in kpc. 

f Corrected HI equivalent line- width at 20% level of intensity (km/s). 

g Corrected total luminosity (in Lbq and Lkq, respectively). 

h Corrected disk central SB (in L^opc -2 and LkqPc~ 2 , respectively). 

1 Corrected disk central (B — K) color. 

J HI mass in M©. 

k References: 1.- Verheijen (1997); Verheijen & Sancisi (2001). 2.- de Jong (1996a, 1996b). 3.- Bell et al. (2000). 4.- de Blok 
et al.(1995, 1996, 2001). 5.- Beijersbergen et al. (1999). 6.- LEDA. 7.- Cox 2000. 8.- Hiromoto et al. (1983). 9.- Gilmore et al. 
(1990). 10.- Kent et al. (1991). 11.- Graham (2002). 
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3.35 


5.51 


N3769 


9.81 


2.73 


9.94 


1.96 


0.57 


10.18 


2.80 


1.91 


0.88 


2.40 


8.36 


N3877 


10.52 


2.92 


9.60 


1.09 


0.11 


10.57 


2.93 


2.55 


0.90 


2.24 


4.12 


N3893 


10.47 


2.89 


10.15 


1.70 


0.32 


10.64 


2.92 


2.71 


0.80 


2.64 


5.32 


N3917 


10.09 


2.36 


9.57 


0.91 


0.23 


10.20 


2.37 


3.21 


0.64 


4.47 


7.42 


N3949 


10.07 


2.89 


9.79 


1.78 


0.35 


10.25 


2.92 


1.70 


0.76 


2.86 


4.40 


N3953 


10.98 


2.96 


9.90 


0.99 


0.08 


11.01 


2.96 


3.95 


0.87 


2.23 


5.83 


N3972 


9.81 


2.62 


9.43 


1.26 


0.30 


9.96 


2.63 


1.88 


0.70 


4.14 


8.39 


U6917 


9.64 


2.01 


9.40 


0.83 


0.37 


9.84 


2.04 


3.13 


0.53 


6.65 


13.62 


U6923 


9.44 


2.21 


9.07 


0.92 


0.30 


9.59 


2.23 


1.85 


0.72 


3.51 


6.50 


N3992 


10.82 


2.80 


10.25 


1.27 


0.21 


10.93 


2.81 


4.56 


0.64 


4.72 


8.98 


U6983 


9.73 


2.06 


9.82 


1.22 


0.55 


10.08 


2.12 


3.78 


0.65 


4.46 


12.68 


N4051 


10.22 


2.79 


9.75 


1.35 


0.26 


10.34 


2.81 


2.38 


0.82 


2.98 


4.47 


N4088 


10.49 


2.85 


10.12 


1.56 


0.30 


10.64 


2.87 


2.99 


0.84 


2.59 


5.32 
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Table 2. Continued. 



Zavala et al.: Luminous and dark matter in disk galaxies 



Name 


LgM s 


T .... b 
Lg^s.O 


LgM g 


L S E 9 


1 


rc 

% 


LgMbar 


Lg£d,o 




(Vd.m/Vt.m) 


M dyn /Mbar 


71 JT IT 

M dyn /L 


N4100 


10.53 


2.92 


10.12 


1.55 





.28 


10.67 


2.93 


2.93 


0.84 


2.73 


5.17 


N4102 


10.55 


3.34 


9.32 


1.31 





.06 


10.58 


3.34 


1.39 


0.91 


1.67 


6.76 


N3718 


10.64 


2.33 


10.01 


0.89 





.19 


10.73 


2.35 


5.42 


0.47 


6.67 


11.35 


N3729 


10.31 


3.11 


9.29 


1.11 





.09 


10.35 


3.11 


1.69 


1.02 


1.91 


5.50 


U6773 


9.01 


1.90 


8.70 


0.70 





.33 


9.19 


1.93 


1.62 


0.68 


3.85 


7.09 


U6818 


9.05 


1.97 


9.13 


1.08 





.55 


9.39 


2.03 


1.94 


0.61 


5.26 


14.72 


N3985 


9.44 


2.66 


8.98 


1.21 





.26 


9.57 


2.67 


1.15 


0.83 


2.95 


3.81 


U7089 


9.28 


1.82 


9.40 


0.92 





.57 


9.65 


1.87 


3.18 


0.72 


4.00 


11.50 


U7094 


8.74 


1.93 


8.92 


1.08 





.60 


9.14 


1.99 


1.55 


1.05 


1.91 


6.32 


N4117 


9.53 


2.22 


9.12 


1.22 





.28 


9.67 


2.26 


1.55 


0.38 


7.80 


22.07 


N4138 


10.40 


3.28 


8.85 


0.94 





.03 


10.41 


3.28 


1.22 


0.84 


1.94 


7.04 


N4218 


8.97 


2.74 


8.63 


1.46 





.32 


9.13 


2.76 


0.61 


0.80 


2.92 


3.42 


N4220 


10.30 


3.16 


9.13 


1.10 





.06 


10.32 


3.17 


1.43 


0.74 


3.15 


5.39 


U128 


10.05 


1.62 


10.16 


0.73 





.56 


10.41 


1.67 


9.57 


0.41 


12.24 


35.99 


F563-V2 


9.54 


2.11 


9.60 


1.18 





.54 


9.87 


2.16 


2.93 


0.87 


2.66 


7.34 


F568-3 


9.85 


1.96 


9.82 


0.95 





.49 


10.14 


2.00 


4.75 


0.59 


5.87 


14.11 


F574-1 


9.86 


1.65 


9.83 


0.82 





.49 


10.15 


1.71 


6.07 


0.58 


4.87 


12.36 


F583-1 


9.11 


1.67 


9.52 


1.43 





.72 


9.67 


1.87 


2.94 


0.55 


5.56 


26.16 


MilkyWay 


10.64 


2.82 


10.09 


1.38 


0. 


.22 


10.75 


2.83 


3.40 


0.71 


3.42 


10.50 


Andromeda 


11.09 


3.07 


10.05 


1.08 


0. 


.08 


11.13 


3.07 


4.26 


0.87 


2.54 


6.83 



a Total stellar mass (disk+bulge) in M Q . 

b Disk stellar central surface density in Mope -2 . 

c Total gas mass in Mq. 

d Disk gas central surface density in Mqpc -2 . 

e Gas fraction (M g /[M g + M s \). 

f Total baryonic mass (M s + M g ) in M . 

g Disk total central surface density in M Q pc~ 2 . 

h Disk (stellar + gaseous) scale radius in kpc. 

1 Disk-to-total maximum velocity ratio. 

J Dynamical-to-baryonic mass ratio at 5/ia- 

k Dynamical mass-to-light ratio in the band B at 5/id(M / L fl u . 



